Dudley Knox Library, NPS 
Monterey, CA 93943 



NAVAL POSTGRADOATE SCHOOL 

Monterey, California 




THESIS 



PRELIMINARY PROPELLER SELECTION USING THE 
WAGENINGEN B-SCREW SERIES AND A 
GENERAL PURPOSE NON-LINEAR OPTIMIZER 

by 

Michael Peter Smith II 
June 1983 

Thesis Advisors: D. Salinas 

G. N. Vanderplaats 

Approved for public release; distribution unlimited. 



1208-/54 



UNCLASSIF IE D. 



SECURITY CLASSIFICATION OF THIS PAGE rWhwt D%tm Enfrmd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE fand Su6mi*) 

Preliminary Propeller Selection Using 
the Wageningen B-Screw Series and a 
General Purpose Non-Linear Optimizer 


S. TYPE OF REPORT 4 PERIOD COVERED 

Master's Thesis; 

June 1983 


6. PERFORMING ORG. REPORT NUMBER 


i. AU^HORfaJ 

Michael Peter Smith II 


8. Contract or grant NUM0ERr«; 


*■ PERFORHINO organization NAME ANO AOORESS 

Naval Postgraduate School 
Monterey, California 93940 


10. program ELEMENT, project, task 
AREA A WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME ANO AOORESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

June 1983 


13. number of paces 

362 


monitoring AGEnCY NAME A AOORESS^// diUmfnt from Controlling OfUcm) 


15. security CLASS, (ot thta report) 

Unclassified 


i5«. declassification/ downgrading 
schedule 



U. OISTRIRUTION STATEMENT (ct thU Rmpott) 

Approved for public release; distribution unlimited. 



17. DISTRIBUTION STATEMENT (ot tt\9 mbmtrmct fitmrmd in Siock 20, i/ dtHmemnt from Rmport) 



It. SUPPLEMENTARY NOTES 



It. KEY WORDS (Continuo on «ld« 1 / md tdontify by block numbmr) 

Wageningen propeller FORTRAN 

propeller strength design optimization 

propeller selection 

ship powering 

20. ABSTRACT (Continuo on rovormm midm if nocmomory mid idmnttfy by block numbmr) 

This thesis presents the use of a general purpose non-linear 
optimization program in the preliminary stage of ship design for 
the selection of a propeller based on methodical series propeller 
test data. The propeller series utilized is the well-known 
Wageningen B-Series. Three (3) "Design Cases", representing the 
thrust, power and matching approaches to powering problems, are 
formulated as FORTRAN subprogram analysis codes for solution by 



OD I JAN*73 1473 COITION OF 1 NOV 6S IS OBSOLCTe 
S/N 0)02- LF-014-6601 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PACE flWi.n Datm Bnfrac- 



1 



UNCLASSIFIED 

security classification of this rage ( Wh*n Data En<ara« 



(20. ABSTRACT Continued) 



the synthesis/optimization program COPES/CONMIN . Designer 
constraints considered are; 

1) diameter limitation 

2) cavitation limit on expanded area ratio using Keller's 
criterion 

3) strength requirement determined by an empirical relation 
and by a method developed by Schoenherr with 
modifications by the author . 

Objective functions considered are maximized open water 
efficiency and minimized propeller blade weight. Optimized 
solutions to specific problems previously presented by other 
authors are obtained and results are compared. 



S' N 0102- LF- 014- 6601 



UNCLASSIFIED 



2 security classification of this RAGEfWhwi D»tm Enfnd ) 



Approved for piablic release; distribution unlimited. 



Preliminary Propeller Selection Using the 
Wageningen B-Screw Series and a 
General Purpose Non-Linear Optimizer 

by 

Michael Peter Smith II 
Lieutenant, United States Naval Reserve 
B.S., St. John Fisher College, 1971 
B.S.E., The University of Michigan, 1976 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1983 



S5^ '75 
C . / 




Dudley Knox Library, NPS 
Monterey, CA 93943 



ABSTRACT 

This thesis presents the use of a general purpose non- 
linear optimization program in the preliminary stage of ship 
design for the selection of a propeller based on methodical 
series propeller test data. The propeller series utilized is 
the well-knovm Wageningen B-Series. Three (3) "Design Cases", 
representing the thrust, power and matching approaches to 
powering problems, are formulated as FORTRAN subprogram 
analysis codes for solution by the synthesis/optimization 
program COPES/CONMIN. Designer constraints considered are: 

1) diameter limitation 

2) cavitation limit on expanded area ratio using 
Keller's criterion 

3) strength requirement determined by an empirical 
relation and by a method developed by Schoenherr 
with modifications by the author. 

Objective functions considered are maximized open water 

efficiency and minimized propeller blade weight. Optimized 

solutions to specific problems previously presented by other 

authors are obtained and results are compared. 
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INTRODUCTION 



I . 

A. BACKGROUND 

The ship design process, in its most rudimentary form, 
has been formulated and tracked by the utilization of the 
classical design spiral (see Figure 1.1) . The design follows 
a convergent helical path past each major milestone "spoke" 
until, after numerous iterative cycles, the final configura- 
tion is "centered" upon. Whether one attempts to segregate 
the principal phases of Preliminary, Advanced and Contract 
Design into separate spirals or combine these phases in 
series along the entire path to the center, it is not long 
before the designer's roughed-out sketches give way to serious 
"number crunching", specifically that of propulsion power 
estimation. 

To estimate the power required to drive the ship through 

the water at its design speed, a decision must first be made 

as to what type of propulsor (i.e., propeller, water jet, 

paddle wheel, etc.) will be used. For the average case, and 

for the discussion that follows, the marine propeller is 

chosen to be the propulsion device. Since 

a ship propeller may be regarded as a transducer 
that converts the rotational power transmitted 
through the shaft into the translational power to 
propel the ship, [Ref. 18: p. 10] 

the selection and design of this device is obviously an 

important factor in the eventual size (weight and power) of 

the ship's propulsion plant. While hydrodynam.icists provide 
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15 



a myriad of theories and techniques to generate a "custom 
built" (i.e., wake adapted) propeller for the ship under 
consideration, their expertise is usually not required in 
the early stages of preliminary design simply because the 
design has not been refined enough beyond gross estimates. 

At this stage, the designer strives to formulate what is 
possible based on previous experience. For preliminary 
power estimation, previous propeller designs (i.e., "stock" 
propellers) and results from methodical series of model 
propellers are analyzed by the designer in order to select 
the "best" available propeller under various conditions posed 
by the problem under consideration. Three examples of typi- 
cal problems encountered in preliminary ship design are: 

1) Given the ship's effective horsepower at a specific 
speed and estimates of hull performance parameters, which 
propeller, as determ.ined by certain principal characteristics, 
will require the least amount of delivered power from the 
propulsion plant? 

2) Given the delivered power from a specific propulsion 
system in terms of torque and revolution rate at the pro- 
peller/shaft interface and estimates of hull performance 
parameters, which propeller will generate the largest effec- 
tive horsepower and speed parameters? 

3) Given a ship's effective horsepower and speed, various 
hull performance parameters, and the propulsion plant's 
delivered power characteristics, which propeller will "match" 
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these requirements at a minimxom amount of weight for a speci- 
fied material? 

(Author's Note: For the sake of brevity, the three 

selection problems just cited will, henceforth, be referred 
to as "Design Case No. 1", "Design Case No. 2" and "Design 
Case No. 3", respectively.) 

For this study, the methodical propeller series method 
is viewed as the designer's choice for preliminary powering 
analysis. One of the m.ost-widely used methodical series data 
on model propellers is the Wageningen B-Screw Series. 

Initially, the results of the series were presented as tabu- 
lations of non-dimensional thrust and torque coefficients 
(K^ and respectively) versus the non-dimensional advance 
ratio (J) for analytical work and as the familiar "Bp-6" and 
"Bu-6" diagrams for design purposes. As "trial & error" 
design methods performed by hand in all engineering disci- 
plines gradually transcended to numerical manipulation by 
the modern digital computer, the necessity for the adaption 
of the Series results to a format suitable for use in computer- 
aided design methods became obvious. This was accomplished 
through multiple regression analysis of the original open- 
water test data of the 120 propeller models in the Series 
and presented in the form of polynomial expressions for "K^" 
and "Kq" [Refs. 1,2]. 

The adaptation of the Wageningen B-Screw Series polynomials 
to various types of propeller selection problems formulated 



17 



for computer solution has been implemented recently by two 
authors. Triantafyllou [Ref. 3] and, of late, Markussen 
[Ref. 4] presented different propeller selection problems and 
proposed different schemes for computer-aided "optimized" 
solutions. In short, specific expressions for the con- 
straints imposed and the objective (optimality condition) 
to be maximized, expressed in terms of a number of design 
variables and parameters, were developed. Then, each system 
of equations was solved by a Newton-Raphson method to give 
a solution set of the design variables which maximized the 
objective and met all constraints. 

Rather than formulating and coding a different optimiza- 
tion scheme each time a propeller selection problem presents 
a different combination and number of design param.eters, 
variables and constraints, a better approach would involve 
formulating the problem (constraints and objective function) 
once in terms of all design parameters and variables and 
utilizing a general purpose optimization schem.e which can 
handle any combination and number of constraints and design 
variables. This alternative certainly allows the designer 
more flexibility in solving his problem. Moreover, it elimin- 
ates repetitive coding and debugging associated with the 
implementation of a computer-sided solution for each particu- 
lar design problem^ 
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B . PROBLEM STATEMENT 

The problem, then, is that the previously cited computer- 
aided "optimized" solutions to the propeller selection problem 
are not broad enough in capability to handle variations in 
the problem formulation. The objective of this thesis is to 
apply an available general purpose optimization com.puter 
code to the solution of various propeller selection problems 
encountered in Preliminary Ship Design in order to enhance 
the flexibility of the selection procedure. 

C. SCOPE 

To achieve the stated objective, the general purpose 
non-linear optimization code CONMIN [Refs. 5,6] together with 
the engineering synthesis code COPES [Ref. 7] (hereafter 
referred to collectively as COPES/CONMIN) is utilized in the 
solution of the three previously cited prelim.inary design 
propeller selection problems. Using the Wageningen B-Screw 
Series propeller characteristics expressed in polynomial 
expressions of various design variables, three "analysis" 
codes, required by COPES/CONMIN, are developed in such a way 
that various combinations of design variables and constraints 
are used, thereby demonstrating the applicability of the 
COPES/CONMIN optimization program in the solution of propeller 
selection problems. 

D. THESIS ORGANIZATION 

The remainder of the thesis is organized in the following 
manner. 
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Chapter II presents a short description of the optimiza- 
tion problem in general terms and a follow-on discussion of 
the COPES/COMMIN optimization program and the mathematical 
techniques employed therein. 

Chapter III introduces definitions and concepts applica- 
ble to the propeller selection problem. A subsequent dis- 
cussion on the Wageningen B-Screw Series is followed by final 
comm.ents on constraints imposed on the propeller selection 
problem.. 

Chapter IV presents the formulation of the propeller 
selection problem as a design optim.ization problem which 
can be solved using COPES/CONMIN. 

Chapter V discusses the background, formulation and 
programming utilized in estimiating a propeller blade's weight 
for subsequent consideration as an objective function. 

Chapter VI reviews the author's modifications to the 
propeller strength analysis developed by Schoenherr [Ref. 8] 
in the early 1960 's for the American Bureau of Shipping. A 
subsequent discussion on the programming details of FORTRAN 
codes, which are utilized for the determination of adequate 
propeller blade strength, com^pletes the chapter. 

Chapter VII reviews the formulation and programming for 
the analysis code which is used in solving propeller selection 
problems represented by Design Case No. 1. Sample solutions 
are presented and compared to those presented previously by 
other authors. 
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Chapters VIII and IX consider Design Case No. 2 and 
Design Case No. 3 selection problems, respectively, in a 
similar fashion to Chapter VII. 

Chapter X, the final chapter, presents the author's 
conclusions and recommendations. 

As a final note, all computer coding presented in this 
thesis is done in FORTRAN IV, the language used by COPES/ 
CONMIN. For the reader's convenience. Appendix A provides 
a cross-reference of the symbols presented throughout the 
thesis to appropriate FORTRAN variable names appearing in 
the author's codes. 
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II. OPTIMIZATION 



A. INTRODUCTION 

The purpose of this chapter is to introduce definitions 
and concepts used in the formulation and solution of the 
general optimization problem. Then, a short discussion on 
the theory and implementation details of COPES/CONMIN is 
presented. 

For further study on the theory and methods of optimiza- 
tion, the reader is directed to the texts by Fox [Ref. 9], 
Fiacco and McCormick [Ref. 10], and Himmelblau [Ref. 11]. 

B. DEFINITIONS 

Before discussing the techniques of optimization and 
their application to engineering problems, some preliminary 
definitions of basic terminology should be stated. Terms 
which have relevant significance are: 

1) Parameters — The numerical quantities for which values 
are assigned to produce a design are called parameters. From 
this, it follows that a design may be specified by a vector 
D containing ”p" components, each of which is associated with 
a param.eter. That is: 



D 

P 
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However, in a design process, the parameters are determined 
by some logical procedure through analysis of some kind. 

Some might take on fixed values to become "preassigned" 
parameters. Interrelationships among other parameters might 
exist so that only some of the parameters are changed when 
one design is compared to another. This consequence leads 
to the definition of "design variable". 

2) Design Variables — The parameters for which values are 
chosen in some fashion to produce a design are called design 
variables. They represent an ordered collection of components 
which is a subset of the design vector D. This subset is 
unique in that its components are "variable", i.e., they may 
take on different values in the design process. Having 
"preassigned" or fixed some of the design's parameters and 
only allowing the remaining "design variables" to change, leads 
to the conclusion that a design is now uniquely specified by 
a vector X containing "n" components (n _< p) , each of which 
is associated with a design variable. That is: 



3) Objective Function — The computable function of all or 
some of the design's preassigned parameters and/or design 
variables, with respect to which the design is to be opti- 
mized, is called the objective function. Single valued in 



X 




( 2 . 2 ) 
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quantitative terms, the objective function's minimum or maxi- 
mum value represents the "best" obtainable or "optimized" 
design. It is expressed as F(D) to show its dependence on 
the design's parameters. But, since a design can be uniquely 
defined by X alone, then clearly F{X) suffices as an expression 
for the objective function. 

4) Constraints — Restrictions on the design which must be 
satisfied in order to produce an acceptable design are 
called constraints. A constraint may be classified as a 
"side" or a "behavior" constraint. A side constraint re- 
stricts or bounds the range of the design for reasons other 
than direct consideration of performance. The side constraint 
on the "i"th design variable may be expressed as: 



lower 

^i 




— 1 



i 



1 , . . . , n 



(2.3) 



A constraint derived from those performance or behavior 
requirements that are explicitly considered is called a be- 
havior constraint. Most often, it appears as a computable 
functional relation involving the design's parameters, both 
preassigned and variable alike. The relation may be an in- 
equality so that the "j"th of "m" inequality constraints 
can be expressed as: 

Gj (D) <0 j=l,...,m (2.4) 
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Alternatively, the relation may be an equality on the "k"th 
of equality constraints expressed as; 



Hj^(D) = 0 k = 1, . . . ,£ (2.5) 

Of noteworthy importance here is that, as before, if 
some of the design's parameters are preassigned, then the 
resulting design is that defined by X which contains only 
the parameters that can be varied in the design process, 
i.e., the design variables. Therefore, constraints imposed 
upon the design may be expressed under one equation as : 



G.(X) < 0 

D 



^ 1,..., m 



H^(X) = 0 



k = !,...,£ 



< X. < 

1 — 1 — 1 



i = 1 , . . . ,n 



( 2 . 6 ) 



A final form for constraints is that of the discrete-valued 
design variable. 

5) Feasible Design — A design in which specified constraints 
are satisfied is called a feasible or "acceptable" design. 

6) Infeasible Design — A design in which constraints are 
violated is called an infeasible or "unacceptable" design. 

C . PROBLEM STATEMENT 

If one presupposes that a range of designs exists within 
a selected design concept, then it follows that different 
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methodologies also exist by which one may choose the param- 
eters which describe the design. One such method is optimi- 
zation where parameters are chosen in a way that the design 
will satisfy all of the limitations and restrictions imposed 
upon it and will be "best" in some sense. In view of the 
foregoing definitions, optimization is then a selection method 
applied to a design problem by which an objective function 
F(D) is minimized to produce an acceptable design which satis- 
fies a certain set of requirements called constraints. 

Formulated mathematically, the general, non-linear, 
constrained optimization problem may be stated under one 
equation as: 



Minimize: F(D) = OBJ 

Subject to: (D) £ 0 

Hj^(D) < 0 

< X. < x“PP®'" 

1 — 1 — 1 



j 

k 

i 



1 

1 

1 



,m 



(2.7) 



Again, as pointed out in the previous section, the design 
may be uniquely defined by just its design variables as 
specified by X when some parameters are preassigned. Thus, 
the general, non-linear, constrained optimization problem 
can now be stated under one equation as: 
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Minimize : 



F(X) 



OBJ 



Subject to: G(X) < 0 j = l,...,m (2.8) 

Hj^(X) = 0 k = 1, . . . ,5. 

< X. < i = 1 n 

1 — 1 — 1 

Solutions methods for this optimization problem, are 
abundant. Those pertaining to the linear and quadratic 
optimization problems involving a few design variables are 
most often presented in graphical or analytic form, although 
numerical schemes are, by no means, a dormant form. Struc- 
tural and thermal problem solutions are most prevalent. How- 
ever, as the optimization problem becomes more complex in 
terms of non-linear relationships among an increasing number 
of design variables and of an increased number of design 
constraints, numerical or mathematical programming techniques 
dominate the solution methods . 

To limit the scope of this discussion, only the numerical 
techniques relevant to COPES/CONMIN will be considered. For 
more background on optimization techniques and applications, 
the reader is directed to a recent paper by Vanderplaats 
[Ref. 12] which presents a concise, but thorough, qualitative 
review of optimization. Although this paper deals exclusively 
with the application of design optimization to structural 
problems, it also contains a very extensive and current list 
of references on general techniques and applications of 
optimization . 
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D. COPES/CONMIN 

As previously stated in Chapter I, COPES/CONMIN is the 
collective acronym for the FORTRAN program utilizing the opti- 
mization code CONMIN and the synthesis code COPES. COPES 
stands for control Program for Engineering Synthesis; CONMIN 
is an acronym for CONstrained function MINimization . 

1. CONMIN 

CONMIN is a FORTRAN program, in subroutine form, which 
solves the general non-linear constrained optimization prob- 
lem as stated; 



Minimize; F(X) = OBJ 
Subject to; (X) £ 0 j 



j^lower 

i 



< X. < x^PP®^ 
— 1—1 



i 



1 , . . * , m 

1, . . . ,n 



(2.9) 



Equation (2.9) applies to the entire statement. Observe 
that equation (2.9) differs from equation (2.8) in that the 
equality constraint set, given by I^(X) =0, is not specified. 
This is because the version of COPES/CONMIN used in this 
study does not consider these types of constraints. However, 
this will not pose any difficulty in solving the propeller 
selection problems previously cited. 

Again, F(X) is the objective function (OBJ). The 
vector X contains the "n" design variables (NDV) . Gj (X) 
are the "m" inequality behavior constraints (NCON) imposed 
on the optimization problem; and x^PP®^ are the 
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respective lower and upper side constraints which bound the 
"design space" over which F(X) and (X) are defined. As 
functional relationships involving X, F(X) and (X) may be 
implicit or explicit, but, in any event, must be continuous 
and have finite numerical values . 

When the inequality condition of equation (2.9) is 
not satisfied, i.e., Gj (X) > 0 for any constraint, the con- 
straint is said to be violated. If the equality condition 
is met, i.e., G^ (X) = 0 for any constraint, the constraint 
is said to be active. And, finally, if the inequality condi- 
tion is satisfied, i.e., G^ (X) < 0, for any constraint, that 
constraint is termed inactive. Any design, defined by X, 
which satisfies the inequalities of equation (2.9) is desig- 
nated as a feasible design. Likewise, any one which vio- 
lates these inequalities is termed an infeasible design. 

The feasible design with the minimum objective function value, 
often referred to as the "minimum feasible design", will, 
therefore, be the optimum design. 

During the optimization process, CONMIN employs the 
Fletcher-Reeves algorithm [Ref. 13] for locally unconstrained 
problems, and Zoutendijk's method of deasible directions [Refs. 
14,15] for locally constrained problems, in a numerical pro- 
cedure which attempts to minimize the objective function, 

F(X) = OBJ, until one or more of the constraints, G^ (X) , 
becomes active. The numerical search procedure begins with 
an initial X vector which may or may not specify a feasible 
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design. Modifications are included in CONMIN so that, if 
the initial design is infeasible, a feasible solution will 
be obtained with minimal increase in F (X) . By iteratively 
updating the design vector X by the following relation: 



^(q+1) ^ ^(q) g(q) 



( 2 . 10 ) 



the optimization process continues by following the con- 
straint boundaries in a direction of search S so that the 
value of F(X) decreases with each iteration q. The scalar 
a* defines the distance of travel in the direction of search 
S. The process terminates when a vector X is found such 
that no further decrease in F(X) can be made. The vector X 
is considered to be optimal and, at least, a local minimum. 

CONMIN can be used alone as a subroutine in any 
FORTRAN program where numerical optimization is desired. 
However, in order to make the optimization process more "user- 
friendly", CONMIN has been coupled to COPES in order to 
simplify its application to various types of problems. 

Further information on CONMIN can be found in previously 
cited references [5] and [6] . 

2. COPES 

COPES is a FORTRAN program that provides automated 
design and trade-off capability to the design engineer. It 
utilizes the optimizer CONMIN to provide the following six 
specific capabilities: 
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.1) simple analysis 

2) optimization 

3) sensitivity analysis 

4) two variable function space analysis 

5) optimum sensitivity 

6) optimization using approximation techniques 

During the execution of COPES, say for optimization, three 
principal tasks are performed; 

1) data management on the design variables and constraints 
through location assignments in a FORTRAN common block called 
GLOB CM. 

2) decision process control on the attainment of an optimal 
design vector X through multiple calls to the optimizer 
until a minimum or maximum value of OBJ is achieved and all 

Gj (X) are satisfied. 

3) evaluation of OBJ and Gj CX) at each X^ and X^"*"^ when 
ICALC = 2 through multiple calls to the user-provided analysis 
subprogram, SUBROUTINE ANALIZ. 

For the application under consideration in this study, only 
the optimization capability will be used. Therefore, further 
elaboration on the other capabilities is not warranted. 

Reference [7] is the user's manual for COPES/CONMIN. 
Details on the mechanics of user implementation are presented 
with subsequent illustration by example. The reader is, 
therefore, encouraged to familiarize himself with the refer- 
ence. However, at this point, it is sufficient to be aware 
of the fact that a user of COPES/CONMIN is required to: 
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1) provide a FORTRAN subroutine called ANALIZ which per- 
forms the input of preassigned parameters, the evaluation of 
the objective function and constraints during the analysis 
phase of the optimization search and the output of the 
results . 

2) provide an assembled deck of control cards required 
by COPES. 

E . CONCLUDING NOTE 

The field of optimization is both extensive and complex 
and, therefore, the foregoing presentation is, by no means, 
complete in every detail. However, it is felt that the pre- 
ceding overview, in conjunction with the cited references, 
covers the necessary prerequisites that will enable the reader 
to follow the application of COPES/CONMIN to the various 
propeller selection problems in the chapters that follow. 



32 



Ill . POWERING, PERFORMANCE AND PROPELLERS 



A. INTRODUCTION 

The purpose of this chapter is to present an overview 
of the terminology and concepts that pertain to ship propul- 
sion, propeller selection and the use of model propeller 
test data. Initially, fundamental definitions used in ship 
powering problems are presented. This is followed by a 
discussion of the "classic" types of propeller design/ 
selection problems encountered by the naval architect and 
marine and naval engineers . Propeller model testing and 
propeller performance characteristics are reviewed next. 

The chapter is completed with a discussion of the VJageningen 
B-Screw Series . 

The goal here is brevity. The reader is, therefore, 
encouraged to investigate the references cited for further 
details . 



B . DEFINITIONS 

Some fundamental terms associated with most propeller 
design/selection problems are; 

1) Effective Horsepower (P^) --power required to tow the 
"bare" hull (without propeller; rudder and appendage allowance 
assumed included) that generates a given resistance (R,p) at 
a given speed (V). It is determined by; 



P 



E 



Rt V 
550 



(3.1) 
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2) Thrust Horsepower (P^) --power delivered to water by 
a propeller developing a thrust force (T) and moving at a 
speed of advance (V^) without the influence of a hull form 
ahead of it. P^j, is determined by: 



P 



T 



T V, 
A 

550 



(3.2) 



3) Delivered Horsepower (Pq) — power delivered by shaft 
to propeller, normally specified at the outboard side of the 
stern tube. Qg is the torque delivered to the propeller; 

Up is the revolution rate of the shaft and, consequently, 
the propeller. Pj^ is determined by: 



P 



D 



2n np 

550 



(3.3) 



4) Shaft Horsepower (P^) — power delivered to the inboard 
side of the stern tube having a transmission efficiency of 
Hg. Pg is determined by 



P 



S 




(3.4) 



5) Brake Horsepower (Pg) — power delivered by the prime 
mover at connection flange to the power train. While Pg is 
normally associated with the prime mover's rated power at 
this connection (BHP) , it can also be specified from the 
power train/propeller side as: 
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(3.5) 



P 



B 




where n n are, respectively, the bearing system and 

D G 

reduction gear transmission efficiencies. 

6) Thrust deduction factor (1-td) — ratio of the tow 
resistance (R,p) to the thrust (T) provided by the propeller. 
It is determined by 



(1-td) 




(3.6) 



7) Wake Factor (Taylor's) (1-wt) — ratio of the speed of 
advance (V^) to the ship's speed (V) . It is given by: 



(1-wt) 




(3.7) 



8) Advance Ratio (J)--a non-dimensional value, associated 
with propeller test data presentation (see figure (3.1), 
given by the following relation: 



V(l-wt) _ 

D D 

P P 



(3.8) 



9) Thrust Coefficient (K^)--a non-dimensional value asso- 
ciated with the thrust force (T) developed by a propeller 
of diameter D which is turning at a rate np and operating 
in a fluid of density p. It is defined by the following 
expression : 
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T 



(3.9) 



K 



T 



2 

n 



10) Torque Coefficient (K^) — a non-dimensional value 
associated with the torque (Qp) absorbed by a propeller of 
diameter D which is turning at a rate np and operating in 
a fluid of density p. It is defined by the following 
expression : 



K 



Q 




(3.10) 



11) Open Water Efficiency (n^) — the ratio of to Pj^ 
for a propeller in open water conditions, with a uniform 
inflow velocity field at a speed of advance V^. It is ex- 
pressed as: 

'’o - np Qg = 2FT^ 



12) Hull Efficiency (p ) — a ratio of work done on the 

n 

ship to that done by the propeller expressed as; 



n 



H 




Rt V 



T V 



A 



(1-td) 

(l-wt) 



(3.12) 



13) Relative Rotative Efficiency (n_) — the ratio of the 
actual, behind-hull efficiency to the open water efficiency. 
The value of Pp does not, in general, depart from the value 
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of 1.0. Most often, ri„ varies between 0.9 5 and 1.0 for twin- 
screw ships and between 1.0 and 1.1 for single screw ships. 
Applicable units for the terms in the expressions above 

are; 

1) horsepower (hp) — P^, P^, Pg , Pj^ and Pg 

2) pounds (Ibf) — R^, T 

3) feet/second (ft/sec) — V, 

4) foot-pounds (ft- Ibf) — Qg , Qp 

5) feet (ft) — D 

6) revolutions/second (rps) — np 

7) revolutions/minute (rpm) — Np = np/60.0 

The quantities T, V^. , D, Q„ and n„ are obtained from the 

It it 

propeller test data results. The quantities R,j, and V are 
specified from the design point on the R-V curve for the hull 
under study. The quantity p is a property of the fluid in 
which the hull and propeller operate. And, finally, n_, 
and Hg are characteristics of the bearing, gear and stern 
tube systems. In preliminary design studies, nominal values, 
based on previous designs, are usually assumed unless, of 
course, these systems have been selected and actual values 
can be specified. 

C . POWERING CONCEPTS 

1 . Basic Relations 

Simply stated, the fundamental powering relationship 
to be solved in ship propulsion and powering problems is: 
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(1-td) 

(1-wt) 



(3.13) 



P 



E 



D 



Utilizing the definitions just presented, equation (3.13) 
can be rewritten as : 



550 



(1-td) 

(1-wt) 



^o 



27T Qg np 
550 



(3.14) 



Rearranging terms of equation (3.14) gives: 
. (l-wt)V _ 

(1-td) 550 ^R ^o * 550 



(3.15) 



And, finally, when substitutions are made, equation (3.15) 
becomes ; 



_ T(l-wt)V 
550 550 



^R ’'c 



277 Qg np 

550 



(3.16) 



Equations (3.14), (3.15), and (3.16) provide the basis for 

different approaches to the solution of a typical powering 
problem. More background and information on the definitions 
and equations presented above may be found in Chapter VI, 
Sections 10-16 in the text by Comstock [Ref. 16] and O'Brien's 
book [Ref. 17] . 

2 . Approaches to the Powering Problem 

From equation (3.16), three types of propeller selec- 
tion problems are discernible. In the first instance, the 
propeller thrust T and the propeller's speed of advance 
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are taken as known quantities. The fact that T is known 
substantiates the "Thrust Approach" nomenclature given to 
this type of selection problem. In the preliminary (or, in 
some circles, conceptual) ship design phase, the specifica- 
tion of T is based upon the requirement imposed by the 
resistance of the ship (R^) at its design speed (V) (or, the 
effective horsepower (P^) at V) and estimates of wt and td 
in the absence of wake surveys and self-propulsion data from 
model tests. Essentially, the thrust delivered by a selected 
propeller must provide, at least, the thrust required for the 
ship hull under study. The objective in the "Thrust Approach" 
selection problem is to determine, by logical means, the 
appropriate values of Qg and np when the open water efficiency 
(Hq) is set by the selected propeller and its performance 
characteristics . 

In the second instance, the delivered torque (Qg) 
the propeller shaft speed (Up) are taken to be known. The 
"Power Approach" nomenclature is given to propeller selection 
problems of this type because P^ is known. Here, with the 
shaft and propeller speeds being equal, the torque absorbed 
by the propeller (Qp) must be, at least, equivalent to the 
delivered torque (Qg) . The corresponding objective in the 
"Power Approach" selection problem is to determine, by logical 
means, the expected ship speed (V) (or, the speed of advance 
(Va) ) and the associated thrust (T) that can be developed 
when the open water efficiency (n^) is, again, set by the 
selected propeller and its performance characteristics. 
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The final, and most familiar, types of propeller 
selection problem occurs when T, or V, Qg and np are all 
known. From equation (3.16), the open water efficiency (n^) 
is now established as a requirement to be met. The objective 
is, simply, to select a propeller whose open water efficiency 
(Hq) , developed thrust and absorbed torque are equivalent to 
or "match" the requirements imposed. Obviously, this approach 
on the selection problem has been designated as a "matching 
problem" . 

The reader is directed to the paper by Vassilopoulos 
[Ref. 18] for further information. 

D. PROPELLER PERFORIVIANCE CHARACTERISTICS 

Up until the late 1950 's, much of the knowledge about 
the performance of propellers has been gained from experience 
with models. To study the relationships governing their 
behavior, a model propeller is built and run in a towing 
tank without any hull ahead of it. This is done by running 
the propeller on a long shaft projecting well ahead of a 
narrow, hydrodynamically shaped pod or "propeller boat" which 
contains the driving mechanism and recording apparatus and 
is attached to the towing carriage. The propeller advances 
into undisturbed fluid (usually water of density p and kine- 
matic viscosity v) so that the speed of advance (V^) is known 
and the flow into the "disc" swept by the turning blades is 
uniform. For the model propeller of diameter (D) under test, 
readings of thrust (T) , torque (Qg) and shaft revolutions (np) 
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are recorded over a range of values for speed of advance (V^) 
in this "open water" condition. 

Using the laws of similitude, the collected data is 
reduced and scaled appropriately into the familiar functional 
relationships between the advance ratio (J) and the non- 
dimensional coefficients of propeller performance. These 
coefficients or performance characteristics, defined pre- 



viously, are 






1) 


Thrust 


Coefficient 


(Kj) 


2) 


Torque 


Coefficient 


(Kq) 


3) 


Open Water Efficiency (n^) 



Figure (3.1) graphically depicts the relationship between J 
and and Kq derived from test data for a propeller defined 
by a specific expanded area ratio (Aj,/Aq) , pitch-diameter 
ratio (P/D) , number of blades (Z) and thickness-to-chord 
ratio (t/c) . 

Definitions of these terms with graphical illustrations 
pertaining to various aspects of propeller geometry can be 
found in Section 15 of references [16], [17] and in van Manen's 

publication [Ref. 19]. 

More recently, highly analytical theories (lifting line, 
modified lifting line, lifting surface, etc.) for use with 
high-speed digital computers have been formulated and subse- 
quently used in "modeling", in a mathematical sense, the 
propeller and its behavior in the "wake adapted" (or, behind 
hull) condition as well as the "open water" condition. 
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Figure 3.1 Open Water Test Results — B 4-100 Series 
Propeller 
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Additional benefits derived from this approach to propeller 
performance analysis include: 

1) determination of blade section profiles along the 
propeller's blade radius (R) to achieve uniform lift and 
internal stress distributions; 

2) computation of "off-design" performance characteris- 
tics in all quadrants; 

3) subsequent determination of hull surface forces, bearing 
loads and spindle torques induced by the propeller; 

4) prediction of steady and unsteady stress distributions 
in the propeller blade using the finite element method on 
the blade of the propeller under study. 

Obviously, this approach to propeller performance analysis 
serves to: 

1) eliminate the time-consioming and expensive model 
construction and testing of propellers in tow tanks and 
cavitation tunnels; 

2) eliminate the "scaling" discrepancies which inhibit 
the reliability of design charts and model propeller data; 

3) eliminate those design charts altogether. 

As in the case with model experiments, however, the ultimate 
objective remains the same, i.e., establishing the performance 
characteristics of the propeller in terms of K^, and as 
functions of J. Having these relationships enables the ship 
designer to proceed in solving the power equation (equation 
(3.13)) through any of the approaches previously discussed. 
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E. THE WAGENINGEN B-SCREW SERIES 



1 . Background 

The model test data of the Wageningen B-Screw Series 
have been selected for use in the powering problems to be 
solved utilizing COPES/CONMIN. The choice was driven by the 
following considerations: 

1) the Series is widely known and, despite its growing 
obsolescence, is still used in preliminary ship design 
studies . 

2) the availability of previous investigations [Refs. 3,4] 
which utilized the series, for comparative analysis of optimi- 
zation results. 

3) the applicability of the polynomial expressions for 
and Kq to computer-aided analysis. 

The Series tests were conducted from 1940 through 1960 and, 
therefore, represent propeller designs (principally naval 
and merchant applications) and design philosophy of that era. 

Specifically, the Series consists of 120 model 
propellers. As is customary in methodical or systematic 
model propeller series testing, the number of blades (Z) , 
expanded area ratio (A^/A^) and pitch-diameter ratio (P/D) are 
varied systematically, while the blade outline, the profile 
of the blade's cross section along the blade radius, blade 
cross section maximum thickness (t) , blade section chord 
length (c) , diameter (D) and propeller hub-to-diameter ratio 
(d/D) were kept constant for given values of A^/A^ and Z. 
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Table (I) summarizes the variations in Z and A„/A^ for each 

L u 

set of model propellers having pitch-diameter ratios (P/D) 
of 0.4, 0.6, 0.8, 1.0, 1.2 and 1.4. 

TABLE I 

Summary of the Wageningen B-Screw Series 




2 . Series Results 

The test results of the Wageningen B-Screw Series were 
originally presented in the form of Bp-6, Bu-5 and K^, K^, 
and -J diagrams. As stated in Chapter I, multiple regression 
analysis was performaed (again, [Refs. 1,2]) on the results 
to produce the polynomial expressions for and Kg. The 
open water efficiency (n^) / as a function of J, follows from 
equation (3.11). The correction for "scale effects" was 
achieved by using Lerb's method of equivalent profiles 
[Ref. 20]. Although Triantafyllou ' s thesis [Ref. 21] sug- 
gests an improved method for scale correction, the results 
of Reference [2] will be used in this study. 
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In propeller selection problems which use this 



Series, values for and Kq are defined as; 



K, 



T 



f^(J,P/D,Ag/AQ,Z,t*/c^^5j^ 



(3.17) 



Kq = f2(J,P/D,Ag/AQ,Z,t*/c^.75j^ 



To compute K,j, and Kq, the following equations are used: 



K, 



T 



+ AK^ 



(3.18) 



K, 



■Q 




The polynomial expressions found in Tables (5) and (6) of 
Reference [2] are then used to evaluate to components K^, 

K^, AK^ and AKq. 

Table (5) in Reference [2] lists the coefficients used 
in the polynomial expressions for and at an equivalent 

g 

Reynolds number (Rn of 2 x 10 . It is defined as: 



Rn 




(3.19) 



.75R 



V 



where ; 



V 



2 

kinematic viscosity of the fluid (ft /sec) ; 



c = blade section chord length at 3/4 propeller 

radius (.75R) in feet (ft). 
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To account for "other effects", coefficients AK^ and AK^ are 
introduced. Table (6) in Reference [2] lists the coeffi- 
cients used in the polynomial expressions for these coeffi- 
cients. "Other effects" include the operation of the 
propeller at an equivalent Reynolds Number different from 
2 X 10 . Also, variations in other parameters which define 
the propeller's geometry, specifically t/c values different 
from the ones fixed by the Wageningen propellers, are taken 
into account by corrections to the equivalent Reynolds number. 
By keeping the blade section's chord length (c) at the value 
of the Wageningen propeller, a change in a blade section's 
maximum thickness from the standard one defined by the Series 
(t ) to one preferred in the selection (t* ) produces a 

new equivalent Reynolds number (Rn*^^^^) given by; 



Rn 



75R 



= exp 



4.6052 + 






(Jin Rn ^^j^-4.6052) 



(3.20) 



where : 



Rn 



* 

.75R 



75R 



Rn 



. 75R 



(t/c) 



.75R 



the new equivalent Reynolds number; 

new equivalent t/c at 3/4 propeller 
radius; 

the Reynolds mamber computed by equation 
(3.19) ; 

standard equivalent t/c at 3/4 propeller 
radius for Wageningen propellers. 
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For the Wageningen B-Screw Series, the standard equivalent 
t/c is given by: 



75R 



(0.0185 - 0.00125Z)Z 
2.073 



(3.21) 



Further details on blade section geometry will be addressed 
in Chapter V. Reference [2] contains background and other 
information on the equations above. 

3. Limitations on Series Data 



propeller selection problem, the following restrictions 
apply to the Series data: 

1) Number of Propeller Blades (Z) — The Series considers 
only propellers with numbers of blades as shown in Table 
( I) . Therefore , 



However, the two bladed propeller, i.e., Z = 2, is not very 
common in conventional merchant and naval ship designs and, 
therefore, is not included in this study. 

2) Equivalent Reynolds Number--The Series data, as pub- 
lished, is valid only in the range of equivalent Reynolds 
nxambers given by: 



In utilizing the Wageningen B-Screw Series in any 



Z 



2 , 3 , 4 , 5 , 6 or 7 



(3.22) 



2 X 10 



6 



< Rn 



< 2 X 10 



9 



(3.23) 



.75R 
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If the equivalent t/c is varied from the standard equivalent 



value (t/c) then the new equivalent Reynolds number 

(Rn*^^j^) , which results from this variation, must lie within 
the same limits. That is. 



2 X 10^ £ Rn*^^j^ £ 2 X 10' 



(3.24) 



These limits are appropriate for full-size propellers. For 
example, given the following; 



a) 


wt 


b) 


V ^ 


c) 


Np 


d) 


D ^ 


e) 


G 


f) 


V ■ 



= 4.0 (ft) 



the value for Rn is equal to 3.619 xlO^. 

3) Pitch-diameter Ratio (P/D) — The series data considers 
only pitch diameter ratios in the range given by: 



0.4 < P/D < 1.4 (3.25) 

4) Advance Ratio (J) — An inspection of the Series results 
in graphical format shows that J varies over a range given 
by; 



0 < J < 1.6 



(3.26) 
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5) Expanded Area Ratio (A^/A^) — Using Table (I) , A^/A^ 
varies over certain ranges depending on Z. This is stated 
as : 



0.35 < A^/Aq £ 0.8 Z = 3 (3.27) 

0.40 < A^/Aq <1.0 Z=4 (3.28) 

0.45 < A„/A^ <1.05 Z = 5 (3.29) 

— iii o — 

0.50 < Ajj/Aq <0.8 Z = 6 (3.30) 

0.55 < A^/Aq <0.85 Z = 7 (3.31) 

6) Hub diameter- to-Propeller Diameter Ratio (d/D)--From 



Table 37, Section 17 of Reference [16], the Series data 
requires that: 



d/D = 0.18 Z = 3,7 (3.32) 

d/D = 0.167 Z = 4,5,6 (3.33) 



F . SUMMARY 

From the preceding discussions, the following observa- 
tions can be made: 

1) the "Design Cases", defined in Chapter I, are examples 
of the powering equation solution approaches. That is. Design 
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Case No. 1 constitutes a "Thrust Approach" problem; Design 
Case No. 2, a "Power Approach" one; Design Case No. 3 , a 
"Matching" problem. 

2) equations (3.17) and (3.8) imply that an optimization 
solution to the "Design Cases" will involve P^, V, wt, D, 

P/D, A^/Aq, (tVc) Qg and np as possible design variables. 

3) when viewed from the concepts on optimization presented 
in Chapter II, equations (3.25) through (3.31) constitute 
side constraints to an optimized solution of a propeller 
selection problem which uses the Wageningen B-Screw Series. 
Having noted these points, the propeller selection problem 
can now be formulated as a general, non-linear, constrained 
optimization problem. 
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IV. PROPELLER SELECTION — AN OPTIMIZATION PROBLEM 



A. INTRODUCTION 

The purpose of this chapter is to present the formula- 
tion of the propeller selection problem as an optimization 
problem that can be solved using COPES/CONMIN. Three usual 
restrictions considered by the designer in any propeller 
selection analysis are stated as constraints. Then, the 
components of D and X are assembled based on requirements 
from previously cited relationships. The restrictions 
considered by the designer and the limitations imposed by 
the use of the Wageningen B-Screw Series are presented in 
inequality constraint format. A formal statement of the 
propeller selection problem as an optimization problem is 
followed by a review of the GLOBCM common block format and 
the basic subprograms used in all three versions of SUBROUTINE 
ANALIZ that pertain to each Design Case. 

The FORTRAN subprogram listings are found in Appendix 

B. Comment cards have been used extensively in the coding 
development to assist the reader. 

B. DESIGNER'S CONSIDERATIONS 

1 . Propeller Size 

The first restriction on the selection of any 
propeller is size. That is, the propeller race in the stern 
of the hull under consideration will only accommodate a 
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As a con- 



propeller of some given maximum diameter . 

straint on a selected propeller of diameter D, this may be 
written as ; 



D 



< D. . 

— lim 



(4.1) 



or, alternatively, as: 

Gg(X) = ^ 1 1 0 (4.2) 

*^lim 

2 . Cavitation 

Another item of importance in propeller selection is 
the cavitation phenomenon. When a propeller of given diameter 
D and expanded area ratio Aj,/Aq is operating to produce a 
thrust T, the formation and subsequent collapse of water vapor 
bubbles on the blade surface, i.e., cavitation, is likely to 
occur if the localized surface pressures, usually on the 
"back" side of the blade, drop below the pressure at which 
the fluid would boil (P^^tvap^ surrounding environment. 

Avoidance of cavitation can be reasonably assured by selecting 
a propeller having certain geometric characteristics . A good 
empirical relationship that establishes these characteristics 
for propellers typified by the Wageningen B-Screw Series is 
the Keller Cavitation criterion [Ref. 2: p. 259]. It speci- 
fies the minim^am required expanded area ratio (A^/Aq 
avoid cavitation and is given by: 
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(4.3) 



(1.3 + 0.3Z) 



(A /A ) = 

' 0 min (P +. + p*acgh „ - p , 

'^atm ^ ^ cl ^watvap 






where: 



Z 

T 

D 

^atm 

^watvap 

P 



acg 



h 

b 



cil 



number of blades; 
developed thrust (Ibf) ; 
propeller diameter (ft) ; 
atmospheric pressure (psia) ; 

fluid vaporization pressure (psia) 

fluid density (Ibf sec^/ft^) 

32.174 (ft/sec^) 

depth to shaft centerline (ft) 

constant: 0.1 for Z = 2, 0.2 for Z = 1. 



As a constraint on the propeller selection, this requirement 
is writen as: 



^^E^'^O^min 






(4.4) 



or 






(A-/A.)„, 

£ ITl.1. XI. 



- 1 



(4.5) 



3 . Strength 

The final designer's consideration (for this study), 
included in the selection of a propeller, is that of strength. 
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Given the propeller's material (promat), selected from Table 

(II) , and the loadings (T and Qg) imposed, it is important to 

ensure that the blade's cross sections have proper dimensions 

(in an ideal sense, maximum blade section thickness (t*) and 

chord length (c) ) to ensure adequate strength. Since the 

use of the B-Screw Series requires that the chord length 

(c) vary as a prescribed function of D, Z and A_/A , as given 

in Table 1 of Reference [2] , the adequacy for strength can 

be determined by an appropriately selected value for blade 

section maximum thickness-to-chord ratio (t*/c) alone. So, 

if t* . is the established minimum blade section maximum 
min 

thickness, then the strength requirement follows from the 
constraint given by: 



t* . 
min 



/c 



(t*/c) 



min 



< t*/c 



(4.6) 



The fact that blade section maximum thickness for the B-Screw 
Series varies linearly with the propeller radius (R) allows 
the strength constraint (equation (4.6)) to be evaluated at 
one section along the radius. This point is chosen to be at 
the 3/4 radius (.75R). Therefore, equation (4.6) becomes; 



(t*/c) 



min — 



(t*/c) 



75R 



(4.7) 



or 



G^^(X) 



(tV0^75 



min 



- 1 



(4.8) 
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Reference [2] suggests the following empirical relation for > 

the minimum required equivalent blade section maximum 

thickness -to-chord ratio (t*/c) : 

. /bR min 



Z /O. 0028+0. 21 



(t*/c) 



(2375-1125P/D)P. 



2 2 

4.123N„D^(S ) 

^ 127788 



min 



2.073 A^/Aq 



(4.9) 



where : 

D = propeller diameter (ft) ; 

Pjj = delivered power (hp) ; 

Np = propeller revolution rate (rpm) ; 

= propeller material allowable stress (psi) ; 
P/D = pitch-diameter ratio. 



However, in Chapter VI of this thesis, an algorithm which 
employs the Schoenherr formulation [Ref. 8] with some modi- 
fications, is presented as an alternative to equation (4.9). 



C. THE DESIGN VECTOR 

In view of the preceding presentations on optimization 
and powering, the design vector D can be assembled for the 
general propeller selection problem utilizing the B-Screw 
Series. This vector is composed of preassigned parameters 
relating to environmental conditions, hull characteristics 
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and the propeller which are required for various equations 
and the design variables. 

1 . Parameters 

a. Environmental 

These parameters pertain primarily to the fluid 
conditions in which the propeller operates and to the atmos- 
phere. Required for various calculations, they are: 

1) fluid temperature (‘’F) — Temp 

2) fluid density (Ibf sec^/ft*^) — p 

2 

3) fluid viscosity (ft /sec) — v 

4) fluid vaporization pressure (psia) — Pv/atvap 

5) atmospheric pressure (psia) — P^^m 

b. Hull Characteristics 

These parameters pertain to certain details pre- 
scribed for the hull under study in the powering analysis. 
They are: 

1) wake fraction — wt 

2) thrust deduction — td 

3) relative rotative efficiency — rio 

4) number of propellers--noscrw 

5) shaft centerline depth (ft) — h . 

O /C 

6) propeller diameter limit (ft) — 

c. Propeller 

These parameters are specified in view of their 
discrete-valued nature. They are: 



1) number of blades — Z 

2) material — promat 
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Table (II) lists materials and properties considered in this 
study. These values are taken from Table (35), Section 15 
of Reference [16] . 



TABLE II 

Material Identifier Reference 



promat 


Material Allowable Stress--Sc 

(psi) 


Density 
(lb f/in 


1 


Cast Iron 


3600—3950 


.260 


2 


Cast Steel 


5915—6265 


.289 


3 


Type 2 Bronze 


7200-7585 


.305 


4 


Type 4 Ni-Al 
Bronze 


8910—9430 


.278 


5 


Stainless Steel 


5400 — 5500 


.283 



2 . Design Variables 

In view of equations (3.13) through (3.17), the 
design variables common to all selection approaches are: 



1) 




2) 


V 


3) 


D 


4) 


P/D 


5) 


V^O 


6) 


(t*/c) 


7) 




8) 


Qs 



The vectors D and X are shown schematically in figure (4.1) . 
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Temp 



D 




P 

V 



*^watvap 

^atm 

wt 

td 

"r 

noscrw 

h n 
cZ 

*^lim 



promat 

V 

D 

P/D 



(t*/c) 



N, 



Figure 4.1 



Design Vectors D and X 



59 



D . CONSTRAINTS 



Besides the constraints imposed by equations (4.2) , 

(4.5) and (4.8), equations (3.25) through (3.31) are rearranged 
to the format of constraints in equation (2.9) . They are 
listed as follows: 

1) Equivalent Reynolds Number — Equation (3.25) becomes: 



1 



.75R 

X 10 ^ 



1000 



(4.10) 



Two constraints are derived; 



Rn 



G3(X) 



= 1 - 



75R 



2 X 10 



6 - 



< 0 



Rn 



G4(X) 



.75R 



2 X 10 



- 1000 < 0 



(4.11) 



(4.12) 



2) Expanded Area Ratio — Equations (3.27) through (3.31) 
become; 



<V''o' lower <"> i V^'o i <V''o> upper 

(4.13) 



Two constraints are derived; 

' 'V*o>lower'^> ' V*0 i ° 

Gg(X) = - (Aj,/Aq)^pp^^(Z) < 0 (4.15) 
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3) Advance Ratio — Equation (3.22) becomes: 



0 < < 1 (4.16) 

— 1.0 — 

Two constraints are derived: 

Gi(X) = < 0 (4.17) 

G2(X) = 3^ - 1 < 0 (4.18) 

4) Equivalent Blade Section Maximum Thickness- to-Chord 
Ratio — Using equation (3.19), boundaries on the range of (t*/c) 
are defined by: 

■^(t/c)^75R £ (t*/c)^^5j^ £ 4(t/c)^.75j^ (4.19) 

Two constraints are derived: 

G7(X) = .75R “ .75R - ^ (4.20) 

Gg(X) = (t*/c) ^. 753 ^ - 4(t/c) < 0 (4.21) 

E. OBJECTIVE FUNCTIONS 

Upon consideration of equation (3.13), Design Case No. 1 
and Design Case No. 2 require that the open water efficiency 
(Oq) / given by equation (3.11), be maximized . In terminology 
related to optimization, this is stated as: 



.75R 
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i'iWi 

'iiiiifi 




(4.22) 



Design Case No. 3, the "matching" problem, requires that the 
blade weight (bldwt) be minimized. This is stated as; 

OBJ^ = bldwt (4.23) 

F. PROPELLER SELECTION OPTIMIZATION PROBLEM STATEMENT 

As a general, non-linear constrained optimization problem 
to be solved by COPES/CONMIN, the propeller selection problem 
for all Design Cases may be stated as one equation given by: 

Minimize: F(X) = OBJ^^ 2 OBJ^ (4.24) 

Subject to: (X) £ 0 j = 1,...,12 

< X. < x“PP®^ i = 1 8 

1 ~ 1 ~ 1 

The constraint G, ..(X) and the values for and 

12 11 

will be specified according to each Design Case. 

G. CODING FUNDAMENTALS 

1 . GLOBCM Common Block 

The GLOBCM common block, required by COPES/CONMIN, 
is now assembled. Table (III) specifies the assignment loca- 
tions for the FORTRAN variables which define objective 
functions, design variables and constraints. 

2. SUBROUTINE ANAL I Z 

While each Design Case uses a different approach, 
all analyses are very similar. Therefore, each SUBROUTINE 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



TABLE III 



Global Common (GLOBCM) Catalog 



FORTRAN DEFINITION 

Name 



ETAO 


^O 




WEIGHT 


bldwt 




AEDVAO 


Ar/Ao 




DIA 


D 




N 


N„ = 60 • n 






P 


P 


PE 


^E 




PDIVD 


P/D 




QS 


Qs 




TC75R 


(t*/c) 




V 


V (ft/sec) 




RJCHL 


Gj^ (X) — eqn 


(4.17) 


RJCNU 


(X) — eqn 


(4.18) 


R75RCL 


G^ (X) — eqn 


(4.11) 


R75RCU 


G^ (X)— eqn 


(4.12) 


AEAOCL 


G^ (X) — eqn 


(4.14) 


AEAOCU 


Gg (X) — eqn 


(4.15) 


TC75CL 


G^ (X) — eqn 


(4.20) 


TC75CU 


Gg (X) — eqn 


(4.21) 


POWBAL 


G, - (X) — eqn 


(7.10 




or 


(9.4) 


DIACNU 


Gg (X) — eqn 


(4.2) 


AEAOCV 


Gio (X) — eqn 


(4.5) 


TCSTRS 


Gii (X) — eqn 


(4.8) 


RJ 


J 





or (8.11) 
(9.6) 
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ANALIZ shares a coinmon structure and other coinmon subroutines 



which perforin calculations required in all cases. Appendices 
C, F and I contain, respectively, the source listings of 
SUBROUTINE ANALIZ for Design Case No. 1, Design Case No. 2 
and Design Case No. 3. 

a. Structure 

The structure common to all cases follows 

accordingly: 

1) all initialization of environmental, hull and propeller 
parameters is accomplished in the input section (ICALC = 1) . 

2) evaluation of and K^, all constraints and appropriate 

objective functions or bldwt) are accomplished in the 

execution section (ICALC = 2) . 

3) output of results for each optimization problem is 
accomplished in the output section (ICALC = 3) . 

b. Basic Subprograms 

The following FORTRAN subprograms are used in 
all three SUBROUTINE ANALIZ codes: 

1) SUBROUTINE CH75RA — calculates the equivalent blade 

section chord length (c for the propeller using Table 1 

[Ref. 2, p. 252] . 

2) SUBROUTINE REY75R — calculates the equivalent Reynolds 
number (Rn*^^^^) using equations (3.19) and (3.20). 

3) SUBROUTINE COEFSA — calculates the thrust and torque 
coefficients (K^ and Kq) through sequential calls to SUBROU- 
TINE CALCKT and SUBROUTINE CALCKQ. The polynomial expressions 
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(Tables (5) and (6), [Ref. 2]) for these coefficients are 
contained in SUBROUTINE CALCKT and SUBROUTINE CALCKQ respec- 
tively. 

4) SUBROUTINE OPWEFF — calculates the open water efficiency 
(Hq) using equation (3.11). 

5) SUBROUTINE JCNA — calculates the constraints on the 
advance ratio (J) given by equations (4.17) and (4.18). 

6) SUBROUTINE REYCNA — calculates the equivalent Reynolds 
number constraints given by equations (4.11) and (4.12). 

7) SUBROUTINE EXTCCN — calculates the constraints on ex- 

panded area ratio (A^/A^) and equivalent blade section maxi- 
mum thickness-to-chord ratio (t*/c given by equations 

(4.14), (4.15), (4.20) and (4.21). 

8) SUBROUTINE DICNUA — calculates the constraint on the 

propeller diameter (D) given by equation (4.2) using the 
hull parameter on maximum diameter . 

9) SUBROUTINE CAVCNA — calculates the constraint for 
cavitation given by equation (4.5) using equation (4.3). 

10) SUBROUTINE STRCNA — calculates the constraint for 
strength given by equation (4.8) using equation (4.9). 

H . SUMMARY 

The propeller selection problem has now been formulated 
as a constrained optimization problem which can be solved 
by COPES/CONMIN. Two items remain for discussion before 
proceeding to specify the final details pertaining to each 
SUBROUTINE ANALIZ code and to present numerical examples. 
These items are; 
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1) the theory and coding relating to the computation of 
the propeller's blade weight for the evaluation of the objec- 
tive function in Design Case No. 3 (OBJ^) . 

2) the theory and coding relating to the computation of 

the minimum required equivalent blade section maximxam thickness- 

to-chord ratio (t*/c -.g. . ) for use in the alternative 

. /D min 

evaluation of the strength constraint given by equation 
(4.8) . 
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V. PROPELLER BLADE WEIGHT — AN OBJECTIVE FUNCTION 



A. INTRODUCTION 

In this chapter, the method for the computation of the 
propeller's blade weight (bldwt) , the objective function 
OBJ, is examined. First, a brief overview on the steps in 
the computational procedure is presented. The FORTRAN 
subprogram SUBROUTINE WGTCAL developed from the algorithm is 
then described. Again, Appendix B contains all subprogram 
listings . 

B. THEORY AND PROCEDURE 

Given the material of a propeller, the calculation of 
the weight of one blade involves nothing more than a volume 
calculation, a relatively routine task performed by most naval 
architects/marine engineers. Analagous to the determination 
of the underwater volume of a ship's hull, the calculation 
is an integration of blade section profiles' cross-sectional 
areas over the propeller radius (R) . 

1. Limits of Integration 

Figure (5.1) depicts a side elevation view of a 
blade and hub, parallel to the propeller shaft axis. The 
cross-hatched area indicates the trace of the volume to be 
calculated. In view of equations (3.32) and (3.33), limits 
of integration are from r= .167Rtor=R for Z = 4,5,6 and 
r = .18R to r = R for Z = 3,7. For convenience, a non- 
dimensional variable "x" will be defined as: 
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1.0 T 



.75 - 



.50 - 



.25 ~ 

.20 

.18 = 
.167 



x=r/R 




Figure 5.1 Propeller Blade & Hub — Side View 
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X 



r 

R 



(5.1) 



Limits are now expressed as x = .167 or .18 to x = 1.0. 



Quite obviously, R = (D/2.0). 



2. Blade Section Profile 



Figure (5.2) depicts an expanded cylindrical blade 



section in profile view at a given r or x. For the Wageningen 
B-Screw Series, the profile is defined, geometrically, by 
a succession of vertical ordinates which specify points along 
the blade section's profile on the "face" (y^) and on the 
back (yj^) with respect to the pitch reference line. At any 
r = xR, vertical ordinates for "aft" (P < 0) and "fwd" 

(P > 0) portions of the blade section are determined by: 




(5.2) 




(5.3) 



and 



yff = v^(t* - t*^) P > 0 



(5.4) 




(5.5) 



where ; 




tabulated values depending on x and P 
(see Tables (2) and (3), Ref. [2]); 
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GENERATOR LINE 
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Figure 5.2 Expanded Cylindrical Blade Section — Profile View 



blade section leading edge thickness (ft) ; 



* 




t* = blade section trailing edge thickness (ft) . 



Units for Y£^/ Yff Yj^£ feet (ft) . 

study, a reasonable assumption is made in that: 



For this 



t 



* 

£e 





t* 



(5.6) 



3 . Blade Section Cross-Sectional Area 

The cross-sectional area at each x = r/R (A(x) ) is 
determined by: 



A(x) 



9 

I 



■1 

1=1 1=1 



10 

I 



A 



fi 



(5.7) 



where 



h . + h . 

A = AP (_ai+i — Si) 
ai ai 2 



(5.8) 



h.. = Yi 



ai 



bai 



- Y 



fai 



(5.9) 



h . 



ai+1 



= Y 



bai+1 



- Y 



fai+1 



(5.10) 



Expressions for A^^/ h^^ and follow in similar fashion. 

Values for y^ . and y, . are determined at 9 points 
^ fai -^bai ^ 

along the "aft" portion of a given blade section's chord (c) : 
values of Y£££ Y£,££ ^.re determined at 10 points along the 
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"fwd" portion. The values for and are fractional 

values of the blade section's chord length (c) at radius 

r = xR as determined from Tables (2) and (3) in Reference 

2 

[2] . The units for A(x) are square feet (ft ) . Units for 

h •/ h ..w h-., h^..,, c, AP . and AP^. are feet (ft), 
ai ai+1' fi' fi+1' ' ai fi 

4 . Volume Integration 

The blade volume (bldvol) is finally determined by 
using Simpson's Rule for integration of A(x) along the non- 
dimensional radius x using appropriate limits. 

5 . Blade Weight 

Once the blade volume (bldvol) is calculated, the 
weight (bldwt) is determined by: - 

bldwt = bldvol • wd .1728 (5.11) 



where: 

bldvol = volume of one blade (ft^) ; 

3 

wd = material weight density (Ibf/in ) . 

Weight Density (wd) depends on blade material (promat) . 

Table (II) lists appropriate values. 

C. CODING 

SUBROUTINE WGTCAL is the main subprogram for the blade 
weight calculation. It, in turn, calls the following FORTRAN 
subprogram for various calculations : 
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1) SUBROUTINE TDIST — generates, at specified radius 
values, a distribution of blade section maximum thicknesses 
(t*) . 

2) SUBROUTINE BLDPRP — generates, at specified radius values 

(i.e., r = .167R or . 18R, . 2R, . 3R, . 4R, . . . , . 9R, 1 . OR) , 

various "blade section properties", one of which is a blade 
section's cross-sectional area given by equation (5.7). 

Other properties which are determined (for later use in 
direct stress computations) include blade section chord 
lengths and centroids and "critical point" locations as 
defined in Chapter VI . 

3) SUBROUTINE BLDVOL — performs a Simpson's Rule integra- 
tion for the propeller blade volume (bldvol) using blade 
section cross-sectional areas generated in SUBROUTINE BLDPRP. 

The blade weight (bldwt) is computed as a final step in the 
main subprogram SUBROUTINE WGTCAL. 

Examination of the codes in Appendix B reveals extensive 
use of common blocks for passing data from one subprogram to 
another. Comment cards provide a full definition of all 
common blocks as well as a description of the task being 
performed at various points in a given subprogram. 

D . SUMMARY 

The coding developed for this study is, admittedly, not 
very compact and efficient. However, the intention has been 
to write all codes with sufficient documentation in order to 
facilitate the reader's understanding of the algorithms employed 
as well as to make the author's debugging work easier. 
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VI. THICKNESS -TO-CHORD RATIO — A DESIGN CONSTRAINT 



A. INTRODUCTION 

The purpose of this chapter is to examine the development 
of an algorithm that will be used to determine the minimum 
required equivalent blade section maximum thickness-to-chord 
ratio used in equation (4.8). The formulation is based upon 
the method developed by Dr. Karl E. Schoenherr [Ref. 8] in 
1963. After a review of the past and present methods employed 
in propeller strength analysis is conducted, a description of 
the Schoenherr model and a list of the assumptions used with 
that model is presented. Then, a brief restatement of his 
model's equations which are used in the algorithm is followed 
by a derivation of the author's modifications to the Schoenherr 
method. The chapter is completed by conducting a review of 
the theory and coding employed by the algorithm. 

The principal reference which is cited throughout this 
chapter is, again. Reference [8]. The reader is encouraged 
to review this reference for further details. 

B. PROPELLER STRENGTH ANALYSIS— A HISTORICAL REVIEW 

Marine propeller blades present a special class of struc- 
tural problem. That problem lies in the difficulty of des- 
cribing a blade design in simple mathematical terms for 
subsequent analysis through various means. Until the "finite 
element method era", analytical methods, including the one 
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[Ref. 8] adapted for this study, relied heavily on practical 
experience of the propeller designer and semi-theoretical 
considerations. Analysis by these methods provided a cri- 
terion of stress rather than actual computation of stresses. 
These methods for predicting blade stresses were developed by 
using "beam" theory or "shell" theory. 

The use of elementary beam theory in propeller strength 
analysis was first adopted by Taylor [Ref. 22] . He treated 
a blade as a cantilever beam attached to the propeller hub 
and loaded by thrust and torque forces distributed linearly 
over the propeller radius. His approach is often deemed a 
"modified beam theory" because he chose to calculate the 
direct stresses using the moment of inertia properties of 
expanded cylindrical blade sections with neutral axis parallel 
to the nose-tail (pitch-reference) line or chord of that 
expanded section. Reasonable estimates of stresses along the 
blade surface were achieved for the unraked, unskewed and 
narrow-bladed propellers of his time. 

As propellers "modernized" and became skewed and wider 
with increasing rake (mostly aft), modifications, improvements 
and alternatives to Taylor's theory were developed. Prin- 
cipally, modifications by Rosingh [Ref. 23] and Hancock 
[Ref. 24] proposed using moment of inertia properties of a 
blade section that was normal to the generating line of the 
axially projected blade outline. Roinson [Ref. 25] later 
improved Taylor's theory for application to wide-bladed 
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propellers. Morgan [Ref. 26] provided an improved method 
for calculating the geometric properties of "modern" airfoil- 
shaped blade sections. Aernoldus and Keyser's [Ref. 27] 
"quasi-static" modeling of the propeller blade allowed 
for additional consideration to stresses induced by centri- 
fugal loading of the raked and skewed blade. The beam 
theory approaches to propeller blade stress analysis culminated, 
for all practical purposes, with Schoenherr's work [Ref. 8] 
in 1963. 

Alternatives to Taylor's beam theory approach, prior to 
1963, consisted of the application of "shell" theory to the 
propeller blade strength problem. This approach was first 
proposed by Conn [Ref. 28] and siibsequently formulated by 
Cohen [Ref. 29] who modeled the blade as a helicoidal shell 
with variable thickness and infinite width. "Shell Theory" 
was utilized again in experimental studies by Connolly 
[Ref. 30] who, like his predecessors, was also forced into 
making an assumption about the behavior of the displacements 
of the blade sections (i.e., constant displacements normal 
to the constant pitch blade at each fractional radius dis- 
tance from the hub) beyond usual assumptions of shell theory. 
Essentially, his experimental results on one specific pro- 
peller contradicted the computational values. Attempts at 
a generalized numerical solution to Connolly's equations 
appeared in 1963 [Ref. 31] and 1964 [Ref. 32]. In 1968, 

Atkinson [Ref. 33], compared Connolly's results with currently 



76 



adopted cantilever beam methods and found, based on the 
inconsistency of results, that it was not possible to, 
recommend one method over the other in the blade strength 
design procedure; another approach was needed. 

Commenting on Atkinson's paper at that time, Sontvedt 
[Ref. 34] pointed out that, in view of these inconsistencies, 
the only approach to blade strength analysis which did not 
require very broad assumptions was the finite element tech- 
nique. Developments in the method at that time were provid- 
ing a new and powerful tool for structural analysis. The 
propeller blade was just another application. Genalis [Ref. 
35] developed codes for the determination of displacements 
and stresses in a blade under hydrodynamic loads using the 
FEM technique and modeling the blade as a shell, a 3-D ele- 
ment mesh of tetrahedrons and rectangular prisms and, finally, 
a composite of shell and 3-D elements. As an aside, a 
finite "difference" solution to Connolly's analytical equa- 
tions was proposed in 1972 [Ref. 36]. In 1973, Atkinson 
[Ref. 37] reported the application of both hydrodynamic and 
centrifugal loads to a blade modeled by a thin-shell triangu- 
lar mesh and a thick-shell parabolic and cubic curved element 
mesh. The results of the triangular element were considered 
unsatisfactory. Another use of the thin-shell triangular 
element was reported by Sontvedt [Ref. 34] in 1974 using the 
SESAM-69 code [Ref. 38] . 

The need to model the blade correctly near the hub, 
where root stresses are usually critical, necessitated the 
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consideration of 3-D elements in lieu of thick/thin-shell 
elements [Refs. 39,40]. The use of the 4- and 10-noded 
tetrahedral element (i.e., TET-4, TET-10 respectively) to 
construct a blade mesh was conducted by Seek [Refs. 41,42]. 

He observed that improved accuracy of stress values, achieved 
by the use of these meshes, were overshadowed somewhat by 
the extensive storage capacity required for each analysis. 
Another "natural" improvement, from a geometrical standpoint, 
appeared in 1978. A general 3-D curved isoparametric element 
was incorporated in a computer code [Ref. 43] developed by 
Ma based on his previous formulation work [Ref. 44] . 

The finite element approach will continue to grow in use 
in propeller blade strength analysis with each successive 
improvement made to the basic elements which are used in the 
mesh generation of the blade. But, until the computer storage 
problem is resolved to the point where one mesh generation 
and subsequent stress analysis of one particular blade be- 
comes a minor processing task, the basic analytical techniques 
will continue to be a meri table "check" [Ref. 45] in the 
preliminary (or conceptual) phase of propeller selection/ 
design*. In this context, the method formulated by Schoenherr 
and his colleagues twenty years ago is considered for adoption 
in determining a minimvim required blade section maximiam 
thickness-to-chord ratio. 
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C . SCHOENHERR ' S METHOD 



1 . Backgroiind 

Schoenherr's method is applicable to preliminary 
(or, conceptual) propeller selection problems because it 
employs an assumed thrust and torque force loading distribu- 
tion for the propeller blade. This assumption is made at 
this stage of the ship's design because the exact wake 
velocity distribution at the propeller race is generally 
not known. 

Also, his method is applicable to propellers repre- 
sented by the Wageningen B-Screw Series because the blades of 
these propellers meet Schoenherr's criteria for the blade 
types covered by his formulation. Specifically, B-Screw 
Series blades have: 

1) a small constant rake angle of 15° over the entire 
blade radius; 

2) a constant pitch distribution over the propeller 
radius with the exception of the Z = 4 propeller whose pitch 
is slightly reduced near the hub; 

3) mild skew; 

4) linear distribution of blade section maximum thickness 
over the radius of the blade; 

5) aerfoil profile qualities where the nose-tail line 

or the chord of the blade section is approximately parallel 
to the pitch reference line. 
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2 . The Blade Model 



Schoenherr models the propeller blade as a cantilever 
beam with unsymmetrical and variable area cross sections 
subjected to loading distributions of hydrodynamic and cen- 
trifugal forces. The following additional assumptions apply: 

1) Flexure theory applies. This subsequently implies the 
following; a) plane cross sections remain plane under load, 
b) Hooke's Law is valid, c) the blade material is homogeneous 
and isotropic, d) fibers are free to extend and contract 
independently of adjacent fibers, and e) stresses at a point 
arising from various forces superimpose. 

2) Shearing stresses and their effects are neglected. 

Only the direct stresses on a strength section are taken 
into account. 

3) The strength sections are taken to be the expanded 
cylindrical blade sections at various radial locations . 

4) The neutral axes of a strength section are straight 
lines passing through the centroid of the expanded cylindrical 
blade section and are parallel and normal, respectively, to 
the pitch reference line, and therefore, the chord, at each 
blade section. 

5) Bending Moments are applied in two planes which are 
mutually perpendicular to each other. One plane is normal 
to the pitch-reference line (and chord) of the strength 
section; the other is parallel. 

6) The angle between the principal axes of inertia and 
the neutral axes is zero . 
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Using this model and assumptions, Schoenherr applies 

the following formula for the evaluation of the direct fiber 

stress ( [a ] ) in a blade section at a radius r = r : 
o o 



to]. 



[M]„^ U [M]„ [F^]^ 

no o ^ Zo o CO 



no 



■£o 



A(x^) 

o 



( 6 . 1 ) 



where : 



[M] 



no ' 



[M] 



Zo 






u 



o' 



w 



o 



A(Xo) 



I 



£o 



I 

no 



X 

o 



resultant bending moments in planes 
normal ( [M] ) and parallel ( [M] ^ to 

the strengtn^section * s chord at r^= r 
(ft-lbf); ° 

centrifugal force acting normal to the 
plane of the strength section at r = r 
and resulting from the centrifugal 
acceleration of the remaining blade 
element mass above that strength section 
(Ibf) ; 

coordinates of a point on the strength 
section's periphery with respect to 
that section's neutral axes system (Z-n 
system) (ft) ; 

strength section's cross-sectional area 
(ft-^) ; 

moment of inertia of the strength 
section with respect to the "JJ," axis 
(ft^) ; 

moment of inertia of the strength section 
with respect to the "n" axis (ft^) ; 

non-dimensional radius given by 

X = r /R. 
o o' 



Since equation (6.1) indicates that the direct fiber stress 
is greatest at points on the periphery of the strength sec- 
tion, Schoenherr selects to examine four "critical points" 
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on the periphery where the fiber stress is likely to be a 
maximxm. These points are designated as 0' ©' © and 

on Figure (6.1) and are specified by coordinates (ul,wl) , 
(u2/w2), (u3,w3) and (u4,w4) respectively in the "£-n" 

reference system. 

The values for [M] _ and [M] „ are determined by 

no £o ^ 

the following relation; 



[M] 

[M] 



no 



£o 



+ 

P no 



'“pile 



[M 

[M 



cb^ no 
cb ^ £o 



( 6 . 2 ) 



(6.3) 



where ; 



[Mp] 



no 



£o 



[M 



cb 



no 



Wcb'jo 



total bending moment due to hydrodynamic 
loading acting in a plane normal to a 
strength section’s chord at r = r 
(ftlbf); ° 

total bending moment due to hydrodynamic 
loading acting in a plane parallel to 
a strength section's chord at r = r 
(ftlbf) ; ° 

total bending moment due to centrifugal 
loading acting in a plane normal to a 
strength section's chord at r = r 
(ftlbf); ° 

total bending moment due to centrifugal 
loading acting in a plane parallel to 
a strength section's chord at r = r 
(ftlbf). ° 



3 . Bending Moments Due to Hydrodynamic Loading 

The derivations of [M ] and [M„] . follow directly 

p no P £o ^ 

from Part I of Schoenherr's paper [Ref. 8: p. 83-89] and. 
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GENERATOR LINE 




o 

u 

il 

5-1 
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Figure 6.1 Strength Section at 



therefore, only key equations will be restated. References 
to equations which contain no decimal point apply to equations 
as numbered in his paper. 

Thrust and torque force are components of the hydro- 
dynamic "lifting” force acting on a blade. Using an assumed 
non-linear distribution of thrust along the blade radius 
given by equation (2) , Schoenherr derives the following ex- 
pression for the bending moment due to thrust ([M^]^) which 

acts at a blade section located at ratius r = r ; 

o 






TR *2*^0* 






where: 



(6.4) 



T 

R 

Z 



= propeller thrust (Ibf) ; 
= propeller radius (ft) ; 

= no. of blades; 






functions of non-dimensional 
radius x evaluated at x = r /R 
and X, = 0.2 and given 6y equa- 
tions^ (4) and (9). 



For the bending 
acts at a blade section 
derives the following; 



moment due to torque ( [M ] ) which 

^ ^ o 

located at radius r = r^, Schoenherr 



[M ] 
q o 



z * 



(6.5) 



where : 
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propeller torque (ft-lbf) ; 
no. of blades; 



Z 

4>p(x = functions of non-dimensional radius x 

° evaluated at x = r /R and x, = 0.2 

and given by equatiSns (19) Snd (4) . 



Figure (6.2) depicts the component resolution for 
[MpJno ^rid which results when equations (6.4) and 

(6.5) are imposed at a strength section at r = r^ which has 



a pitch angle 3 



The following relations are derived as 



equations (42) and (43) in Schoenherr's paper: 



‘'^p'no = ‘“t’o * ‘“q'o 

4 . Force and Bending Moments Due to Centrifugal Loading 
The derivation and expressions contained in this 
section constitute the author's modifications to the formu- 
lation in Part II of Schoenherr's paper. In Part II, Schoen- 
herr's derivations for [M , ] and [M , ] „ are formulated for 

cb no CD Zo 

computation using a propeller drawing. This follows from 
the fact that his method, which was funded by the American 
Bureau of Shipping, was intended to be used as that classi- 
fication society's "designer's check" on adherence to the 
Bureau's strength criteria from a propeller blueprint. To 
evaluate the direct fiber stresses from equation (6.1) for 
the Wageningen B-Screw Series in accordance with Schoenherr's 
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GENERATOR LINE 
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Figure 6.2 Bending Moments due to Thrust and Torque 



method/ [M , ] and [M , ] „ must be evaluated from expres- 
cb no cb zo 

sions derived from the available information on blade sec- 
tion profiles and other geometric characteristics which 
are contained in Reference [2] and previously used in Chapter 

V. 

Consider Figure (6.3) where the centrifugal force 

[C„] of a blade element above a blade section at x = x = r /R 
F o o o' 

acts in a radial direction from the shaft centerline. Its 

line of action passes through point "N", which is on the same 

cylindrical surface as the blade section at x = x^ = 

and through point "G" , which is the center of gravity of the 

blade element above the blade section at x = x = r /R. 

o o' 

From the figure, the following expression is derived: 

^ sin i;(x^) (6.8) 

[C„] is shifted to point "N" and is decomposed into 

£ O 

components [C„] cos r(x ) and [C„] sin c(x ) . The entire 

r O ^ O r O O 

cylindrical surface in which the blade section at x = x^ = r^/R 
and point "N" lie is now expanded into a flat plane for further 
consideration (see Figure (6.4)). In this configuration, 

[C„] cos r(x ) is normal to this flat plane while [C„] sin c(x ) 

r O O r O O 

lies in this plane. 

Let point ”0" be the location of the blade section's 
neutral axes system (i.e., the £-n system). Then, the forces 
and moments due to the centrifugal reaction of the blade ele- 
ment above this section act at point "0" and are given by: 



87 



PITCH REFERENCE LINE GENEPATOR LINE 




li 

X 

o 







88 



Figure 6.3 Center of gravity "G" and, intersection point 



GENERATOR LINE 
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Figuce 6.4 Internal Loading on a Blade Section at 



where : 



[F ] 
CO 






[M , ] 

cb o 



[M ] 

cw o 



direct force, due to centrifugal action, 
acting on the blade section located at 

^ = ""o = 

shear force, due to centrifugal action, 

acting on the blade section located 

at X = = r /R; 

o o 

bending moment at point "0" imposed by 
[F^]^ acting through point ”N"; 

torsional moment at point "0" imposed by 
[D^]^ acting through point "N". 



Since Schoenherr's method does not consider shear forces and 

their effects, [D_] and [M ] will not be considered in 

c o cw o 

this modification. However, [F ] and [M , ] must now be 

CO cb o 

computed for each blade section along the propeller's radius 

in order to account for their contributions to equations 

(6.2), (6.3) and, finally, in equation (6.1). 

The computation is derived as follows. Consider 

Figure (6.4). Again, [F ] acts through point "N" and is 

normal (outward) to the plane of the figure. [M , ] is now 

cb o 

resolved into components of the l-n axes system as follows; 



(6.13) 



no 



[Cj.]oCOsC(Xo){Po 



sin B + q cos B^ + Y } 
O o o CO 



5 (x ) {q sin B„-p„ cos B„ + X } (6.14) 

cb Zo F O O O O O O CO 



where ; 



B = pitch angle of the blade section at 

° X = X = r /R; 

o o' 

X = distance of the blade sections centroid from 

CO 

the generator line (ft) ; 

Y = distance of the blade section's centroid from 

CO 

the pitch reference line (ft) ; 

q = distance to point "N" from point "P" 

parallel to the shaft axis at x = x = r /R 
(ft); ' ° ° 

p = distance to point "N" from point "P" 

° perpendicular to the shaft axis at 

X = X = r /R (ft) . 
o o 



The quantity B^ is found by the relation: 



tan B^ = i • (P/D)^ 

O TT O 



(6.15) 



For the Wageningen B-Screw Series, (P/D)^ is a constant along R 
except for propellers with Z = 4 . 

The quantity [C„] is computed from the relation; 

r O 



F o 



= 1728 



wd V 
o 

acg 



( 2Trnp) 



R 



(x_) 

G o 



(6.16) 



where; 
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V 

o 



acg 

wd 



(x ) 

G o 



R 



volume of the blade element above the blade 

section at x = x = r /R (ft^) ; 

o o 

32.174 ft/sec^; 

material weight density (Ibf/in^) ; 

propeller revolution rate (rps) ; 

non-dimensional radial position of "G" 
for the blade element above the shaft 
axis ; 

propeller radius (ft) . 



At this point, only five quantities remain to be determined 
for the evaluation of the expressions of equations (6.10), 

(6.13) and (6.14). They are: 

1) cos c(x^) 

Po 

''o 

These quantities are determined by integration over the 

blade element above the blade section, located at x = x = r /R, 

o o 

from x=x tox=1.0. 
o 

The values for p^ and q^ will vary with the location 
of "G” (from Figure (6.3)) which depends on x. Consider a 
radially thin slice of the blade element above the blade 
section, located at x = x^ = (see Figure (6.5)) . This 

thin "slice" is located at a non-dimensional distance x from 
the shaft centerline where x^ < x < 1.0. Figure (6.6) depicts 
this section expanded onto a plane. Let "g" be the centroid 
of that "slice". If x^ is the distance from the generator 
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Figure 6.6 Coordinates of centroid ’’G" of any Blade Section 




i 

I 








line to "g" in feet and is the distance from the pitch 
reference line in feet, then, from Figures (6.5) and (6.6) , 
the following relations apply; 

i = X R tan r\ (6.17) 

g 

i -a = i -X sin 3 - y cos 6 (6.18) 

g g g g g ^g g 



t = X cos 3 - y sin 3 (6.19) 

g g g ^g g 

where; 

n = rake angle at x; 

X = distance of "g" from point "P" parallel to 
^ the shaft axis (ft); 

y = distance of "g" from point "P" perpendicular 
^ to the shaft axis (ft) . 

For the Wageningen B-Screw Series, rake angle n is a constant 
15° everywhere along the radius R. 

Now, to compute the volume of the blade element above 
the blade section located at x = x^ = r^/R, the following 
expression is used; 



1 

V = f RA(x) dx (6.20) 

o 

x=x 

o 

To compute the non-dimensional radial position of 
"G" for the blade element above a blade section located at 
X = x^ = r^/Rf the following expression is used; 
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1 

/ A (x) X dx 

(^G>o = 4-^ \i.21) 

/ A(x) dx 

x=x 

o 

To compute the tangential position of "G” for the 

blade element above a blade section located at x = x = r /R, 

o o 

the following expression is used; 



1 

/ A ( X ) t dx 

x=x ^ 

T = (6.22) 

O 1 

/ A(x) dx 

x=x 

o 

And, finally, to compute the axial position of "G" 
for the blade element above a blade section located at x = x 

( 

= the following expression is used; 



1 

/ A (x) ( i -a ) dx 
X ^ ^ 

A^ = — (6.23) 

o 1 

/ A{x) dx 

X 

o 

Using the values just determined for 
the following expressions are used to evaluate and 

^o 



P 



o 



X 

o 



<='g>o 



T 

o 



(6.24) 
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X tan n 
o 



(6.25) 



q = A - 
o 



The expression for cos C follows: 



cos 



X R 
o 



(6.26) 



The formulation is now complete. Equations (6.10), 

(6.13) and (6.14) can now be evaluated for any blade section 

located at x = x^ = r^/R. From here, equations (6.2) and 

(6.3) are evaluated. Finally, using the results from these 

equations and equation (6.9), equation (6.1) can be evaluated 

for the four points, specified by coordinates (ul,wl) , 

(u2,w2) , (u3,w3) and (u4,w4), at any location x = x^ = r^/R. 

From the development discussed in Chapter V, the 

values for A(x) , x , y , I and are readily available. 

g g no Zo ^ 

D. ALGORITHM FOR THE CONSTRAINT 



1 . Theory 

Schoenherr's formulation with the modifications just 
derived can be used in determining the minimum required 
equivalent blade section maximum thickness- to-chord ratio 



((t*/c) 



. 75R min 



) for use in the constraint Gj^q^(X) £ 0 given 



by equation (4.8) . The procedure employed is as follows: 

1) assume an initial value for (t*/c) ..j-t. • using 

equation (3.21); 

2) increase this value by a small amount; 
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3) using (t*/c) °t)tained from step (2) , generate 

a distribution of minimum required blade section thicknesses 
(t^in) for blade sections at specified points along the 
propeller radius, say at r = ,2R, ,3R, .4R, .5R, .6R, .7R, 

. 8R and . 9R; 

4) determine all blade section properties to include: 

a) cross-sectional area, b) chord length, c) centroid location, 
d) moments of inertia with respect to the principal axes 
system (i.e., Jl-n system), e) coordinate values for the four 
critical points defined in the previous section; 

5) compute the hydrodynamic bending moment components 

[MpJno radius locations just specified; 

6) compute the values of the centrifugal force [F ] and 

c o 

the bending moment components and acting on 

blade sections at radius locations just specified; 

7) calculate the direct fiber stresses at all four 
critical points for all radius locations specified in step 
(3) ; 

8) check the following condition on the calculated fiber 
stress at all four critical points at all specified radius 
locations us ing ; 



[a] < 144 • S (6.27) 
o c 

9) if the maximum allowable stress (S^) for the material 
is exceeded, then return to step (2) and repeat. Otherwise, 
proceed to next step. 
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10) since the minimiojn required equivalent blade section 
maximum thickness -to-chord ratio assumed in step (2) has 
produced blade sections of adequate strength, evaluate the 
constraint given by equation (4.8). 

2 . Coding Details 

The algorithm just outlined is incorporated into the 
main FORTRAN subprogram SUBROUTINE STRCNK. This subprogram, 
in turn, executes the algorithm through sequential calls to 
other key FORTRAN subprograms. These subprograms are listed 
as follows: 

1) SUBROUTINE TDIST — accomplishes step (3); generates, 

at specified radius values , a distribution of minimum required 

blade section maximum thicknesses (t*. ) for the assumed 

min 

value of (t*/c) 

. /dR min 

2) SUBROUTINE BLDPRP--accomplishes step (4) ; described 
previously in Chapter V; 

3) SUBROUTINE HYDLD — accomplishes step (5); computes the 
hydrodynamic bending moment components , given by equations 
(6.6) and (6.7), at specified radius locations; 

4) SUBROUTINE CNFGLD — accomplishes step (6); computes the 

centrifugal force and bending moments , given by equations 
(6. 9), (6.13) and (6.14) respectively, at specified radius 

locations ; 

5) SUBROUTINE SIGNDS — accomplishes step (7) ; computes 
direct fiber stresses, given by equation (6.1) , for all four 
critical points at every specified radius location. 
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During the remaining steps of SUBROUTINE STRCNK, the 
condition on allowable stress, given by equation (6.27), 
is checked at all critical points of blade sections located 
at specified radius locations (again, r = .2R, .3R, .4R, ,5R, 

.6R, .7R, . 8R and .9R). The final calculation made is 
that for the constraint given by equation (4.8). 

Again, extensive use of common blocks, for passing 
data from one subprogram to another, is apparent upon examina- 
tion of the codes just cited. Comment cards are used 
throughout . 

E . SUMMARY 

The end of this chapter marks the completion of all 
prerequisite background and formulation discussions on the 
application of COPES/CONMIN to propeller selection problems 
involving the Wageningen B-Screw Series. From this point, 
each specific Design Case can now be solved as an optimiza- 
tion problem. 
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VII. DESIGN CASE NO. 1 — PROGRAMMING AND COMPARISONS 



A. INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the "thrust" approach. 
First, the thrust approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solution 
to this problem is presented. Four variations to this pro- 
peller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 
the results from the four variations. 

B. THRUST APPROACH FORMULATION 

1 . Design Vector xT 

As previously pointed out at the conclusion of Chap- 
ter III, Design Case No. 1 constitutes a propeller selection 
problem which is solved by the thrust approach. In this ap- 
proach, the effective horsepower (Pg) and the ship's speed 
(V) are specified by the designer. From the viewpoint of 
optimization, the quantities P^ and V become preassigned 
parameters. This reduces the design vector X (see Figure 4.1) 
to : 
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Having specified and V, all of the design variables, 
as listed in equation (7.1), are not independent. Recalling 
equations (3.3) and (3.13) , the following relationship 
results ; 



P 



E 



(1-td) . 

(1-wt) ^o 



2^Qg Np 

550 60 



Rearranging terms, this equation becomes: 



X, _ (1-wt) 550-60 

^o ^P (1-td) ' T]^ ' 2tt 



(7.2) 



(7.3) 



Considering that the open water efficiency (O^) is_ evaluated 
prior to the computation of (-n ) , or OBJ, -, then both N 

O X ^ Z ir 

and Qg are not independent design variables. One must 
be selected as the independent design variable. Then, 
the other variable becomes dependent on the one just 
selected. 

For this study, Np is selected as the independent 
design variable. This choice will reduce the design vector 
X for propeller selection problems using the thrust approach 
to the following; 

D 

P/D 

<"*/=>. 75R 
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Finally, equation (3.8) implies an alternative defini- 
tion of X as given in equation (7.4) . The design vector for 
Design Case No. 1 propeller selection problems is, therefore, 
defined as: 



! D J 

P/D / 

Ae/Aq > (7.5) 

(t*/c) ( 

J I 

2 . Powering Constraint 

Having determined the design vector XI, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint 0 ^ 2 ^^) men- 
tioned in Chapter IV. 

Simply stated, the selected propeller, as defined by 
XI, must develop enough thrust (T) so that the powering require- 
ment, specified by Pg and V, is met. Using equation (3.2) , 
the thrust developed by the propeller can be specified in terms 
of thrust horsepower (P^) as: 



^^T^ dev 



T V(l-wt) 
550 



(7.6) 



From equation (3.9), it follows that: 



^^T^dev 



P Hp D^ 
550 



V(l-wt) 



(7.7) 
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Using equation (3.12), the developed thrust horsepower can 
be defined in terms of developed effective horsepower given 
by: 



<Vdev = • ‘"’T'dev '’-S' 

The restriction imposed by the thrust approach method, 
where P^, and V are specified, can now be stated as: 






dev 



(7.9) 



Rearranging equation (7.9), the constraint follows: 

dev 

G,-(X1) = 1 ^ Q (7.10) 

^E 

With the design vector xT and defined, the 

propeller selection problem represented by Design Case No. 1 
can be stated under one equation as : 



Maximize: F(X1) = OBJ^ 2 

Subject to: G^ (jTT) £ 0 j = 1,...,12 (7.11) 

< XI. < i = 1 5 

1 — 1 — 1 

C. PREVIOUS SOLUTIONS 

Triantafyllou [Refs. 3,21] considered a propeller selection 
problem represented by Design Case No. 1. In his example 
problem, the following parameters were specified; 
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1) V = 1.139 xio"^ (mVsec) = 1.22613 x 10"^ (ft^sec) 

2) wt = .22 

3) td = .19 

4) = 1.025 

5) noscrw = 1 

6) Z = 5 

7) Pg = 18153 (hp) 

8) V = 24 (knots) 

9) D = 22.0 (ft) 

10) Aj./Aq = .85 

The hull under study in his example had the following dimen- 



sions 


2 




1) 


Length = 710 (ft) 


2) 


Draft = 30 


(ft) 


3) 


Beam = 100 


(ft) 



For his analysis, the design vector contained two varia- 
bles and was specified as: 




Using an iterative scheme [Ref. 21: p. 79] to solve two equa- 
tions in two unknowns, he maximized the open water efficiency 
(tIq) to obtain the following results: 

P/D = 1.1651 
Np = 104 (rpm) 
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25544 (hp) 



n = .6676 

o 

For future comparisons, equations (3.8) and (3.3) give: 

J = .8286 

Qg = 1290000.0 (ft-lbf) 

Triantafyllou ' s results are summarized in Table (IV). 

D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
(7.11), is now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 

The first and second variations attempt to reproduce the 
solution given by Triantafyllou. The design vector XT 
(NDV =2) is used in both cases. One variation uses SUBROUTINE 
STRCNA to evaluate the constraint given by equation 

(4.8). The other uses SUBROUTINE STRCNK to determine • 

The remaining two variations will solve the propeller 
selection problem using the design vector XI (NDV = 5) defined 
in equation (7.5). Again, one variation uses SUBROUTINE 
STRCNA; the other, SUBROUTINE STRCNK. 

In all variations, the following parameters are used: 

1) Temp = 59 (°F) 

2) p = 1.9384 (Ibf-sec^/ft"^) 

3) V = 1.2285 10 (ft^/sec) 
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“> Pwatvap = •2'>7 

Patm = 

6) wt = .22 

7) td = .19 

8) Dj^ = 1.025 

9) noscrw = 1 

10) h = 19.0 (ft) 

11) D, . = 22.0 (ft) 

lim 

12) Z = 5 

13) promat = 5 (stainless steel; see Table (II)) 

14) = 18153 (hp) 

hi 

15) V = 24.0 (knots) 

All of the above are initialized in the input phase (ICALC = 1) 
of each SUBROUTINE ANALIZ pertaining to each variation. 

The constraint evaluated by SUBROU- 

TINE BLPOWl which appears in the execution section of each 
SUBROUTINE ANALIZ. 

1 . Variation 1 

a. Prograinining Details 

Since this variation uses the design vector XT, 
the following design variables of XI become parameters and 
are specified in the input section of SUBROUTINE ANALIZ (ICALC = 1) 
as ; 

1) D = 22.0 (ft) 

2) A^/Aq = .85 

3) (t*/c) = .0348 (from equation 3.21). 
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For constraints, the following are used; 

Gj (XT) <_ 0 j = 1, ... ,8, 12 

Only nine of twelve constraints are evaluated (NCON = 9) . 
Obviously, some of the twelve constraints are redundant since 
D, A^/Aq, and (t*/c) have been specified. 

The upper (XtVP^®^) and lower (XT^®'^®^) limits on 
the design variables J and P/D are set to be: 

.01 < J £ 1.6 

.4 £ P/D <1.4 

These upper and lower limits are specified in the COPES 
control card deck on card image F under respective fields 
VUB and VLB. The initial value for each design variable 
(XT^) is also assigned on card image F under the field labeled 
X. The first list of card images in Appendix D lists all of 
the COPES control cards used for this variation and variation 
2. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7] . 

An examination of SUBROUTINE ANALIZ for this varia- 
tion, found in Appendix C, shows the calling statement made 
to SUBROUTINE STRCNA. 
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b. Results 



The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed first in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 

2 . Variation 2 

a. Programming Details 

Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix C, shows 
the calling statement made to SUBROUTINE STRCNK. 

b. Results 

The output from the optimization/ analysis , per- 
formed by COPES/CONMIN, is listed second in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 

3 . Variation 3 

a. Programming Details 

This variation uses the design vector 30?. For 
constraints, the following are used: 

Gj (3a) < 0 j = 1, ... ,12 

All twelve constraints are evaluated (NCON = 12) . 

The upper (Xl^^^®^) and lower limits 

on the design variables D, P/D, A^/A^, (t*/c) - 75 j^» snd J 

are set as : 
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1.0 1 D < 50.0 (ft) 

.4 £ P/D 1 1.4 
. 2 ^ ^ 

.003 < (t*/c) < .50 

.01 ;< J _< 1.6 

These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (Xl^) 
is also assigned on card image F under the field labeled X. 

The second list of card images in Appendix D lists all of 
the COPES control cards used for this variation and variation 
4. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7] . 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix C, shows the calling statement 
made to SUBROUTINE STRCNA. 
b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN , is listed third in Appendix E. Re- 
sults for this variation of thepropeller selection problem 
are tabulated in Table (IV) . 
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4 . Variation 4 



a. Programming Details 

Everything discussed above for the third varia- 
tion applies here with one exception. An examination of 
SUBROUTINE ANALIZ for the fourth variation, found in Appendix 
C, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed last in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV) . 

E. DISCUSSION 

Overall, the results achieved in all variations compare 
reasonably well to the solution obtained by Triantafyllou. 
However, the following points can be made. 

Variations 1 and 2 give the same results. This was ex- 
pected in view of the fact that, even though constraints 
Gg(XT) through G^j^(XT) were evaluated, these constraints were 
not considered in the optimization search conducted by CONMIN. 

In variations 3 and 4, the diameter (D) was driven to 
the limit ♦ This bears out a fundamental rule in pro- 

peller design, i.e., the larger the propeller diameter (D) , 
the greater the open water efficiency (n^) . 

The minimum required equivalent blade section maximum 

thic]aiess-to-chord ratio ( (t*/c) . ) , computed in 

' . 75R min ' 
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variation 3, is substantially smaller than the one computed 
for variation 4. As pointed out in Reference [2], the empiri- 
cal relation, expressed by equation (4.9) and derived from 
equation (70.x) [Ref. 46: p. 620] , does not take into account 
the effects of centrifugal loading. These effects include, 
specifically, the direct stresses imposed by the inertia load 
of the blade and the bending moments which result from rake 
and skew of the blade. Therefore, the algorithm developed in 
Chapter VI should, and does, produce a larger value for 
(t*/c) 75j^. 

A final observation on the results concerns the values 
of the open water efficiency. The ''optimum” open water effi- 
ciency (Hq) achieved by Triantafyllou is lower than those 
achieved in variations 1 and 2. A possible reason for this 
might be the neglection of the term ”dRe/dJ" in Triantafyllou ' s 
formulation of the analytical expressions [Ref. 21: p. 71] 
that he used in his analysis. The difference in the open 
water efficiencies siobsequently accounts for the differences 
in the propeller revolution rate (Np) and the delivered 
torque (Qg) when the relation in equation (3.3) is considered. 
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TABLE IV 



Design Case No. 1 — Results 



GROUP 


ITEM 


TRLANTA- 

FYLLOU 


1 


2 


3 


4 






18153.0 


18153.0 


18153.0 


18153.0 


18153.0 


Given 


V 


24.0 


24.0 


24.0 


24.0 


24.0 


Design 


D 


22.0 


22.0 


22.0 






Variable 

Speci- 


V\> 


.85 


.85 


.85 






fied 


.75R 




.0348 


.0348 








D 








21.9991 


21.9659 


Design 


P/D 


1.1651 


1.0036 


1.0036 


.9981 


1.0071 


Variables 










.8205 


.8149 












.0330 


.0642 




J 


.8286 


.7371 


.7371 


.7343 


.7394 


Maximize 


'>0 


.6676 


.7091 


.7091 


.7109 


.7066 




^lim 




(22.0) 


(22.0) 


22.0 


22.0 


Restric- 

tions 


^S/\)min 
(tVc) ^75pi;nin 


(0.5258) 

(0.21266) 


(0.5258) 

(0.50761) 


.5269 

.021998 


.5267 

1 .053259 






104 


116.9 


116.9 


117.4 


116.7 


Other 


Qs 


1290000 


1080451. 


1080451 


1064574 


1084003 




^D 


25544.0 


24051.3 


24051.3 


24030.3 


24094.8 
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VIII. DESIGN CASE NO. 2— PROGRAMMING AND COMPARISONS 



A . INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the "power" approach. 
First, the power approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solu- 
tion to this problem is presented. Four variations to this 
propeller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 
the results from the four variations. 

B. POWER APPROACH FORMULATION 
1 . Design Vector X2 

As previously pointed out at the conclusion of 
Chapter III, Design Case No. 2 constitutes a propeller selec- 
tion problem which is solved by the power approach. In this 
approach, the delivered torque (Qg) and the propeller revolu- 
tion rate (Np) are specified by the designer. From the view- 
point of optimization, the quantities Qg and Np become pre- 
assigned parameters. This reduces the design vector X (see 
Figure (4.1)) to; 

^E 
V 

D 

P/D 
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Having specified Qg and Np , all of the design varia- 
bles, as listed in equation (8.1) , are not independent. Re- 
calling equations (3.3), (3.11) and (3.13), the following 

relationship results; 



(1-td) 

(1-wt) 



'R 



J 

2iir 



2ttQ 

w 



60 



( 8 . 2 ) 



Rearranging terms, this equation becomes; 



11 = (1-td) 

V (1-wt) 



'R 



(1-wt) 
np D 



K, 



550 



N 



P 

60 



(8.3) 



Considering the relations for and in equation (3.17), 
then both and V are not independent design variables. 

One must be selected as independent, while the other becomes 
dependent on the one selected. 

For this study, V is selected as the independent de- 
sign variable. This choice will reduce the design vector X 
for propeller selection problems using the power approach to 
the following; 



X = 



V 

D 

P/D 

(t*/c) 



(8.4) 



75R 



Finally, equation (3.8) implies an alternative defini- 
tion of X as given in equation (8.4) . The design vector for 
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Design Case No. 2 propeller selection problems is, therefore, 
defined as; 



X2 



V 

J 

P/D 

(t*/c) 




(8.5) 



2 . Powering Constraint 

Having determined the design vector X2, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint men- 

tioned in Chapter IV. 

Simply stated, the selected propeller, as defined by 
X2, must absorb at least all of the power delivered to it 
(Pj^) which is specified in terms of Qg and Np . Using equation 
(3.3), the power absorbed by the propeller can be specified 
in terms of delivered horsepower (Pq) as: 



^^D^ absorb 



2tt Qp Np 

550 ' ^ 



( 8 . 6 ) 



From equation (3.10), it follows that: 



^^D^ absorb 



*=0 I’ "p 2Tt Nj 

550 ' 60 



(8.7) 



But, equation (3.3) also defines the power delivered to the 
propeller as; 
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2tt N 



550 



P 

60 



( 8 . 8 ) 



The restriction imposed by the power approach method, 
where Q„ and N are specified, can now be stated as: 

^ It 



absorb 



(8.9) 



Rearranging equation (8.9), the constraint Gj^2(^2) follows: 



G^2 (X2) 



f p 1 

1 _ ^ absorb ^ 



( 8 . 10 ) 



Further simplification of equation (8.10) gives 



Gj_2(X2) 



1 - ^ < 
Qs - 



( 8 . 11 ) 



With the design vector X2 and Gj^2^^2) defined, the 
propeller selection problem represented by Design Case No. 2 
can be stated under one equation as; 



Minimize: F(X2) = OBJ^ 2 

Subject to; G ^ (j^) 0 j = 1,...,12 (8.12) 

jj2lower ^ ^2 . < i = 1 5 

1 — 1 — 1 

C. PREVIOUS SOLUTIONS 

Markussen [Ref. 4] considered a propeller selection prob- 
lem represented by Design Case No. 2. In his example problem. 
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the following parameters were specified; 



1) 


Temp =18 (°C) =64.4 (°F) 




2) 


Wvap = 0-0206411 (bars) = 


.29943921 (psia) 


3) 


P ^ = 1.01312856 (bars) = 14 

atm 


.6974 (psia) 


4) 


noscrw = 1 




5) 


h^^ =6.7 (meters) = 21.9827 


(ft) 


6) 


Z = 6 




7) 


Pj^ = 18.9 (Megawatts) = 25344 


.9 (hp) 


8) 


Np = 110 (rpm) 




9) 


= 15.65 (knots) 





For his analysis, the design vector contained three 
variables and was specified as; 

- } 

XM = / P/D V, 

( ) 

A restriction for the minimiam required expanded area ratio 

( (A_/A ) ), given by equation (4.3), was also considered. 

^ u inin 

This imposed a constraint given by equation (4.4). 

Using an iterative scheme [Ref. 4; p. 110] to solve three 
equations in three unknowns, Markussen maximized the open 
water efficiency (ri^) to obtain the following results; 

J = .61095 

P/D = .864380 

Ae/Aq = 36.1861/40.6123 (m^/m^) 
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.891012 



n, 



o 



.654391 



For future comparisons, equations (3.8) and (3.3) give: 



D 



7.19091 , (meters) = 23.593375 (ft) 



Q 



S 



1210130.0 (ft-lbf) 



Markussen's results are summarized in Table (V). 

D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
(8.12), is now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 

The first and second variations attempt to reproduce the 
solution given by Markussen. The design vector XM (NDV = 3) 
is used in both cases. One variation uses SUBROUTINE STRCNA 
to evaluate the constraint 0^2 given by equation (4.8). 

The other uses SUBROUTINE STRCNK to determine 0^2 • 

The remaining two variations will solve the propeller 
selection problem using the design vector X2 (NDV = 5) de- 
fined in equation (8.5). Again, one variation uses SUBROUTINE 
STRCNA; the other, SUBROUTINE STRCNK. 

In all variations, the following parameters are used: 

1) Temp = 64.4 (°F) 



2 ) 



3) V 



p = 1.9892 (Ibf-sec^/ft"^) 

V = 1.1900 x 10“^ (ftVsec) 
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5) 



Patm 14.69 7 (psia) 

6 ) wt = . 2 2 

7) td = .19 

8) n„ = 1.025 

9) noscrw = 1 

10) = 21.9827 (ft) 

11) D, . = 30.0 (ft) 

1 im 

12) Z = 6 

13) promat = 5 (stainless steel; see Table (II)) 

14) Qg = 1210130 (ft-lbf) 

15) Np = 110 (rpm) 

All of the above are initialized in the input phase 
(ICALC = 1) of each SUBROUTINE ANALIZ pertaining to each 
variation. 

The constraint Gj ^2 Is evaluated by 

SUBROUTINE BLPOW2 which appears in the execution section of 
each SUBROUTINE ANALIZ. 

1 . Variation 1 

a. Programming Details 

Since this variation uses the design vector XM, 
the following design variables of 3^ become parameters and 
are specified in the input section of SUBROUTINE ANALIZ 
(ICALC = 1). The ship's speed (V) is specified as: 

V = V^/(l-wt) 

= 15.65/(1 - .22) 

= 20.0641 (Icnots) 
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Markussen elected to use the standard Wageningen blade sec- 
tion maximum thickness-to-chord ratios. Since the equivalent 
t/c is given as a function of Z and Aj,/Aq (see equation (3.21)), 
then (t*/c) is calculated during each analysis (ICALC = 2) 

by the following relation: 

(t*/c)^75R = (t/c)^75R 



For constraints, the following are used: 

Gj (XT) <0 j = 1,...,8, 10, 12 



Only ten of twelve constraints are considered (NCON = 10) . 
Constraints Gg(XM) and G^j^(XM) are redundant since no limit 
on the propeller diameter appears as a parameter in 

Markussen 's formulation and (t*/c) was taken to be the 

Wageningen standard. 

The upper (XmV^P®^) and lower limits 

on the design variables J, P/D and A_/A_ are set to be: 

ii» o 



.01 


< 


J 


< 


1.1 


.4 


< 


P/D 


< 


1.4 


.4 


< 


Ar/Ao 


< 


1.1 



These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VIB. The initial value for each design variable (XM^) is 
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also assigned on card image F under the field labeled X. 

The first list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 2. 
These cards also specify the locations of the design varia- 
bles in the common block GLOBCM (see Table (III) ) as well as 
the locations of the constraints and their boundaries . Further 
details on the COPES control card requirements and the format 
of each card are contained in Reference [7] . 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, 
performed by COPES/CONMIN, is listed first in Appendix H. 
Results for this variation of the propeller selection problem 
are tabulated in Table (V) . 

2 . Variation 2 

a. Programming Details 

Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix F, shows 
the calling statement made to SUBROUTINE STRCNK instead of 
SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed second in Appendix H. 
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Results for this variation of the propeller selection problem 
are tabulated in Table (V) . 



3 . Variation 3 

a. Programming Details 

This variation uses the design vector X2 
constraints, the following are used; 



For 



Gj (X2) £ 0 j = 1, ... ,12 



All twelve constraints are evaluated (NCON = 12) . 

The upper (X2V^^^^) and lower (X2^°^®^) limits 
on the design variables V, P/D, A„/A_, (t*/c) and J are 



set as; 



10.0 


< 


V 


< 


100.0 


.4 


< 


P/D 


< 


1.4 


.4 


< 


^e/^0 


< 


1.1 


.003 


< 


(t*/c) 


< 


.50 


1 — 1 
o 

• 


< 


J 


< 


1 — 1 

• 

1 — 1 



These upper and lower limits are specified in' the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X2^) 
is also assigned on card image F under the field labeled X. 

The second list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 4 . 
These cards also specify the locations of the design variables 
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in the conunon block GLOBCM (see Table (III)) as well as the 
locations of the constraints and their boundaries. Further 
details on the COPES control card requirements and the format 
of each card are contained in Reference [7] . 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN , is listed third in Appendix H. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (V) . 

4 . Variation 4 

a. Programming Details 

Everything discussed above for the third varia- 
tion applies here with one exception. An examination of S’UB- 
ROUTINE ANALIZ for the fourth variation, found in Appendix 
F, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis , per- 
formed by COPES/CONMIN, is listed last in Appendix H. Results 
for this variation of the propeller selection problem are 
tabulated in Table (V) . 

E. DISCUSSION 

The results achieved in variatio ns 1 and 2 compare ex- 
tremely well to the solution obtained by Markussen. As 
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pointed out in the discussion in Chapter VII, variations 1 
and 2 are expected to give the same results for the vector 
3 cm. Obviously, the values obtained for J and P/D, as well 
as those for D, and Rn*^ 2 j^/ very close to the values 

generated in Markussen's example. However, the values for 
A„/A_, are somewhat different. It is interesting to note that 
the value obtained in variations 1 and 2 (and, for that matter, 
variations 3 and 4) is, essentially, the limiting value for 
as given in Table (I) , for Z = 6. Markussen's value 

b U 

for Aj,/Aq (i.e. , .891012) exceeds the limit (i.e., .80) in 
this table. 

As pointed out at the end of Chapter VII, the minimum 
required blade section maximum thickness-to-chord ratio 
( (t*/c) . )) , computed in variations 1 and 3, is sub- 

stantially smaller than the one computed for variations 2 
and 4. Again, the same explanation applies here as well. 

The results of variations 3 and 4 differ somewhat from 
Markussen's results. The reason for this is simply that V^ 

(or V) has not been specified as a parameter. Consequently, 
a higher value for the advance ratio (J) , which corresponds 
to a higher open water efficiency (tIq) / has been found in 
the optimization search. This result can be interpreted 
in the following way. Given; 

1) a six-bladed Wageningen propeller (Z = 6) which is made 
out of stainless steel (promat = 5) ; 

2) a power train delivering 25344.9 (hp) at a rate of 
110 (rpm) ; 
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3) a hull with a wake fraction (wt) equal to .22, a 
thrust deduction (td) equal to .19 and a shaft centerline 
depth (h^^) of 21.98 (ft), 

then, the selected propeller, as defined by X2 , can drive 
this hull at a maximum speed of V when the hull has a maximum 
resistance given by P„. 

ill 
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TABLE V 



GROUP 



Given 



Design 

Variable 

Specified 



Design 

Variables 



Maximize 



Restrictions 



Other 



Design Case No. 2--Results 



ITEM MARKUSSEN VARIATIONS 

12 3 4 



Pd 


25344.9 


25344.9 


25344.9 


25344.9 


25344.9 




1210130 


1210130 


1210130 


1210130 


1210130 




no 


no 


no 


no 


no 




15.65 


15.65 


15.65 






V 




20.0641 


20.0641 






J 


.61095 


.6475 


.6475 


.9927 


.8753 


V 








30.9138 


29.5355 










24.n27 


23.0377 


p/d 


.864380 


.9036 


.9036 


1.1986 


1.0308 


V\) 


.891012 


.8018 


.8018 


.7946 


.7986 


(t*/c) 




.0397 


.0397 


.0499 


.0638 




.654391 


.6660 


.6660 


.7643 


.7330 


^lim 




(30.0) 


(30.0) 


30.0 


30.0 


"^S/^imin 


.574729 


.5070 


.5070 


.4622 


.4266 


(t*/c) 


(.02729) 


(.0647) 


.02706 


.0638 


D 


23.53 


22.25 


22.25 


22.36 


24.23 


Pe 




(14057.3) 


(14057.3) 


19945.5 


20653.7 


^!?5R 


, 478x10”^ 


5x10^ 


5 X 10 


5 X 10 


5x 10 
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IX. DESIGN CASE NO. 3 — PROGRAMMING AND COMPARISONS 



A. INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
a propeller selection problem where "matching" is desired. 
First, the "matching" approach to the propeller selection 
problem is formulated. Then, a review of a previous author's 
solution is presented. One variation to this propeller 
selection problem is solved by COPES/CONMIN. The chapter is 
completed with a presentation and discussion of the results. 

B. "MATCHING" FORMULATION 
1 . Design Vector 

Design Case No. 3, the final powering problem con- 
sidered in this study, constitutes a propeller selection prob- 
lem solved by the "matching" approach. In this approach, 
the hull's effective horsepower (Pg) and speed (V), the 
delivered torque (Qg) and the propeller revolution rate 
(Np) are specified by the designer. This reduces the design 
vector X (see Figure (4.1)) to; 

I D 

P/D 

V^O 

75R 

For this study, the design vector X is reduced further 
by eliminating the propeller diameter (D) as a design variable. 




128 



That is, D will also be specified by the designer so that 
the design vector for Design Case No. 3 is defined as: 

( P/D 

X3 = I A^/A^ 

( (t*/c) 

2 . Powering Constraint ( s) 

Having determined the design vector X3, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint men- 

tioned in Chapter IV as well as an additional constraint. 

In the "matching" problem, the selected propeller, 
as defined by X3 , must satisfy two conditions. First, it 
must develop , as a minimum, the effective horsepower (Pg) 
as imposed by the design specification. Citing the formula- 
tion previously derived in Chapter VII, this condition can 
be stated as : 



(9.2) 



'■e i '^E’dev 



(9.3) 



The constraint follows accordingly as: 



Gj_2(X3) 



^ dev ^ 



0 



(9.4) 



For the second condition, the selected propeller can 
only absorb , as a maximum, the delivered power (P^) as 
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specified by the designer. The formulation is the same as 
that in Chapter VIII except that the inequality signs are 
reversed. The condition is stated as: 



absorb - 



(9.5) 



In defining a constraint another location, say 

location 24, in the GLOBCM block (see Table (III)) would be 
assigned. But, considering the fact that constraint G^ (X) 
will not be used because the propeller diameter (D) is speci- 
fied, there is no reason why G^{X3) cannot be redefined, for 
this Design Case only, as: 

G (X3) = - 1 < 0 (9.6) 

D 

Further simplification of equation (9.6) gives: 

Gq(X3) = - 1 < 0 (9.7) 

9 Qg 

In reality, the constraints just defined should be 
equality constraints. The word "match" does infer equality 
in some sense. However, as previously stated in Chapter II, 
the version of COPES/CONMIN used in this study does not 
directly handle equality constraints. But, since CONMIN 
attempts to minimize constraints in the optimization search, 
it will be assumed that a "match" can be achieved. 
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With the design vector X3 and the constraints 
Gt o (^) G-(3r3) defined, the propeller selection problem 

represented by Design Case No. 3 can now be stated under 
one equation as : 

Minimize; F(X3) = OBJ 3 

Subject to; G^ (X3) ^0 j = 1,...,12 (9.8) 

< X. < X3“PP®'^ i = 1 3 

1 — 1—1 

C. PREVIOUS SOLUTIONS 

The propeller selection problem considered by Vassilopoulos 
[Ref. 18] actually represents a propeller "design" problem 
using the "power" approach. The example problem which he 
elected to solve is taken from that posed by the International 
Towing Tank Conference (ITTC) Propeller Committee. This 
problem is concerned with the determination of propeller 
thrust (T) , diameter (D) and speed of advance (V^) (or ship 
speed (V)) for a single-screw cargo ship where; 

1 ) power available to the propeller (i.e., Pj^) is 
30,000 (hp) 

2) Z = 6 

3) N = 105—110 (rpm) 

4) = 23 (ft) 

5) = 19 (ft) 

The variation of ship speed (V) and of hull effective power 
(Pg) , thrust deduction factor (td) and the wake fraction 
(wt) is also given [Ref. 18; p. 20]. 
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The results from Vassilopoulos ' propeller design exercise 
produced a propeller that is "matched" at the following values 

1) Pg = 21292.6 (hp) 

2) V = 24.24 (knots) 

3) Qg = 1500606.75 (ft-lbf) 

4) Np = 105 (rpm) 

5) Pq = 30000.0 (hp) 

His propeller "design" was based upon the following specified 
parameters : 

1) Temp = 59 (°F) 

2) p = 1.9905 (Ibf-sec^/ft^) 

3) p . = . 247 (psia) 

^watvap ^ 

4) wt = .22 

5) td = .1725 

6) noscrw = 1 

7) h^^ = 19.0 (ft) 

8) Z = 6 

9) promat = 5 (stainless teel, see Table (II)) 

10) D = 22.0 (ft) 

11) Np = 105 (rpm) 

12) Pj^ = 30000.0 (hp) 

Using an optimization scheme incorporated in his MVAPDP 
computer program, Vassilopoulos maximized the open water 
efficiency (n^) and designed a propeller with the following 
characteristics : 

1) J = .852 

2) = .242 
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.0478 



3) Kq = 

4) Hq = .691 

5) A^/Aq = .767 

6) bldwt = 7617.2 (Ibf) 

By utilizing both the lifting line and lifting surface 
methods in his design procedxire, Vassilopoulos ' MVAPDP program 
evolved a "constant stress" propeller blade. Consequently, 
the values for (t*/c) and P/D varied non-linearly along the 
propeller radius (R) . According to Vassilopoulos, this 
resulted in a minimum weight propeller. The values for P/D 
and (t*/c) are listed in Tables (8) and (10) of his paper. 

From these values, (t*/c) is approximately .040. 

While the propeller represented by Vassilopoulos' design- 
is different, in many aspects (rake, skew, blade section 
aerfoil shape, etc.), from the Wageningen B-Screw Series 
propeller, it does represent a minimum weight propeller that 
has been "matched" to specific design values. Appropriate 
results are summarized in Table (VI) . 

D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
(9.8), is now solved by COPES/CONMIN. One solution variation 
is considered. The following parameters are used: 

1) Temp = 59 (,°F) 

2) p = 1.9905 (Ibf-sec^/ft'^) 

3) V = 1.2817x10"^ (ftVsec) 

Pwatvap = *247 (psia) 
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Patm * fPsia) 

6) wt = .22 

7) r]^ = 1.025 

8) noscrw = 1 

9) = 19.0 (ft) 

10) Z = 6 

11) promat = 5 (stainless steel, see Table (II) ) 

12) D = 22.0 (ft) 

Two problems are examined. Problem 1 specifies the following 
additional parameters : 

1) td = .1725 

2) P„ = 21292.6 (hp) 

3) V = 24.24 (knots) 

4) Qg = 1500606.75 (ft-lbf) 

5) Np = 105 (rpm) 

Problem 2 specifies the same parameters as: 

1) td = .171 

2) Pg = 17630.0 (hp) 

3) V = 23.0 (knots) 

4) Qg = 1500606.75 (ft-lbf) 

5) Np = 105 (rpm) 

All of the above are initialized in the input section 
(ICALC = 1) of similar versions of SUBROUTINE ANALIZ. There- 
fore, only one version is included in Appendix I. 

The constraints for Gg {X3) and G^2 evaluated by 
SUBROUTINE BLPOW3 which appears in the execution section of 
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SUBROUTINE ANALIZ. Also, note that SUBROUTINE DICNUA has 
been deleted from the execution section, while SUBROUTINE 
WGTCAL has been added. 

1 . Programming Details 

All twelve constraints are evaluated (NOON = 12) . 

The upper and lower limits on the design 

variables P/D, Aj,/Aq and (t*/c) are set to be; 



.4 


< 


P/D 


< 


1.4 


.4 


< 




< 


1.1 


.003 


< 




< 


.50 



These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X3^) 
is also assigned on card image F under the field labeled X. 
The list of card images in Appendix J lists all of the COPES 
control cards used for both problems . These cards also 
specify the locations of the design variables in the common 
block GLOBCM (see Table (III) ) as well as the locations of 
the constraints and their boundaries. Further details on 
the COPES control card requirements and the format of each 
card are contained in Reference [7] . 

2 . Results 

The outputs from the optimization/analysis, performed 
by COPES/CONMIN, are listed in Appendix K. Results of both 
problems are tabulated in Table (VI) . 
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E. DISCUSSION 



Table (VI) presents the results of problems 1 and 2 
along with relevant information from Vassilopoulos ' "design”. 
Problem 1 attempted to "match" a Wageningen propeller at the 
design point found by Vassilopoulos. The first COPES/CONMIN 
printout in Appendix K indicates that the "match" was achieved 
at Pg equal to 21.168.1 (hp) and P^ equal to 28,150.0 (hp) 

(or, Qg = 1500607 (ft-lbf) and Np = 105 (rpm) ) . These values 
are judged to be close enough to the "Given" values in Table 
(VI) . 

It is apparent that the Wageningen propeller does not 
require all of the 30,000 (hp) of delivered horsepower. The 
propeller characteristics (i.e., J, K,j,, Kq and n^) for problem 
1 compare very well to Vassilopoulos' values. The expanded 
area ratios (Ag/A^) are, also, very similar. Of course, the 
obvious difference is the blade weight (bldwt) . The Wageningen 
propeller blade is over five thousand pounds heavier. Does 
this make sense for a minimum blade weight? 

The answer is yes. 

All one has to do is consider the values of (t*/c) for problem 
1 and Vassilopoulos' design. Vassilopoulos' "constant stress" 
blade was designed to "absorb" stress up to the allowable 
design limit of 5,400 (psi) (for stainless steel) all along 
the entire propeller radius (R) . Table (12) in Reference 
[18] gives further details. The Wageningen propeller blade, 
however, represents an "older" type of blade which was de- 
signed with a linear blade section maximum thickness (t*) 
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distribution. Consequently, it was "overdesigned" for strength 

beyond the 3/10 — 4/10 radius (i.e., . 3R — .4R) and contains 

excess material. A heavier blade, therefore, results. Note, 

also, that the optimizer did not drive the value of (t*/c) 

to the minimum acceptable value, (t*/c) -,ct< 

. / dR min 

The results of problem 2 show the effect on blade weight 
(bldwt) for a Wageningen propeller when the hull's powering 
requirem.ents (i.e., at V) have been reduced. The weight 
reduction of 2000 pounds is significant. The complete re- 
sults are listed in the second COPES/CONMIN printout in 
Appendix K. 
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TABLE VI 





Design 


Case No. 3- 


-Results 




GROUP 


ITEM 


VASSILOPOUIOS 


PROBLEM 








1 


2 






21292.6 


21292.6 


17630.0 


Given 


V 


24.24 


24.24 


23.0 




Qs 


1500607 


1500607 


1500607 






105 


105 


105 


Design 

Variable 


D 


22.0 


22.0 


22.0 


Specified 


P/D 


* 


1.1813 


1.0906 


Design 

Variables 




.767 


.7944 


.7742 






.040 


.0794 


.0681 


Minimize 


bldwt 


— 


12842.6 


10464.7 


Maximize 




.691 


— 


— 




^lim 


— 


— 


— 


Restric- 

tions 


■^V^^min 
(t*/c) _ 


— 


.8515 

.0691 


.7722 

.0681 




J 


.852 


.8290 


.7866 






.242 


.2349 


.2063 


Other 


""q 


.0478 


.0448 


.0372 




% 


— 


.6915 


.6950 




bldvt 


7617.2 


— 


— 



★ 

P/D varies with R 
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X. CONCLUSIONS AND RECOJXIMENDATIONS 



A. CONCLUSIONS 

The general purpose non-linear optimizer/synthesizer 
COPES/CONMIN has been successfully applied to three typical 
preliminary ship design propeller selection problems in which 
the Wageningen B-Screw Series is used. The formulation and 
programming of each required analysis code (i.e., SUBROUTINE 
ANALIZ) have been made as general as possible to allow the 
designer a broad variety of solution options for solving 
propeller selection problems which can be classified under 
any of the three Design Cases that were considered. The 
analysis codes have been "modularized" to the extent that 
methodical series data from other propeller series, which 
are available in the polynomial expression format of the 
B-Screw Series, can be easily adapted for powering analysis 
utilizing design optimization methods. 

Further flexibility in the solution to the propeller 
selection problem has been achieved by using COPES/CONMIN 
as the optimizer/synthesizer. The designer has now been 
afforded the additional capability of specifying the design 
variables, the objective functions and the constraints of his 
choice. By solving propeller selection problems in the way 
presented in this thesis, repetitive problem formulation and 
coding have been eliminated. 

There are other advantages to solving propeller selection 
problems specifically with COPES/CONMIN which have not been 
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directly addressed in this study. As stated in Chapter II, 
COPES/CONMIN is capable of performing optimization analyses, 
sensitivity studies, optimum sensitivity studies and optimi- 
zation using approximation techniques. The designer, there- 
fore, can select and perform any of these options, using the 
same analysis codes which have been presented in this thesis. 

While the utilization of a general purpose non-linear 
optimizer in solving propeller selection problems allows the 
designer greater flexibility in the selection procedure, there 
is one important limitation that should be stressed at this 
point. This concerns the question whether or not the solution 
vector, determined by the optimizer, is a "global" optimum. 

As stated in Chapter II, COPES/CONMIN assures that, if a 
feasible solution vector is found, it is, at least, a "local" 
minimum (or maximum) . This implies that, for two different 
initial design vectors which are specified in the COPES Con- 
trol Card deck on card image F, the same optimum solution 
may not be determined by the optimizer. Both solutions would 
correspond to minimums (or maximums) of the objective function 
and are, therefore, correct. But, does one or the other 
correspond to the minimum (or maximum) of the entire vector 
design space, i.e., the "global" optimum? For the moment, 
at least, there is no definitive answer to the question. 

Despite this uncertainty, progress in the field of design 
optimization continues to be made. Current developments 
[Ref. 47] will soon allow the designer to have a choice in 
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selecting a specific optimization algorithm from a "library" 
of proven optimization programs which employ the latest 
state-of-theart numerical techniques. Again, using one 
analysis code, the designer will be able to generate any 
number of optimized solutions for the problem under study. 

B . RECOMMENDATIONS 

For future consideration, it is recommended that the 
automated design and trade-off capability, provided by a 
general purpose non-linear optimizer/synthesizer such as 
COPES/CONMIN, be applied to the more difficult problem of 
propeller design . 

As pointed out in Chapter I, the use of the Wageningen 
B-Screw Series represents a "selection" procedure rather 
than a "design" process. Today, analytical propeller design 
procedures, utilizing lifting line and lifting surface 
theory, are becoming increasingly popular among propeller 
designers. The propeller design, which results from the 
utilization of these analytical methods, is, unquestionably, 
more efficient than the standard series propeller. However, 
these methods require consideration of many more design varia- 
bles in the design process. This appears to be a natural 
application for the use of a general purpose non-linear 
optimizer/synthesizer. 

Here, an analysis code, much larger than those which 
have been presented in this study, could be developed which 
would incorporate the lifting line/lifting surface theory 
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for the determination of the propeller performance character- 
istics, the local cavitation numbers and also the calculation 
of the pressure distributions over the blade. These pres- 
sure distributions would be utilized in the strength analysis 
of the blade. This analysis would utilize the finite element 
technique on an appropriately generated mesh model of the 
blade. Having defined the steps for this design procedure 
in the analysis code, the propeller designer now "couples" 
his analysis to the optimizer/synthesizer for determination 
of the optimum design. A massive amount of computer storage 
would certainly be required, but this concept is feasible 
and, in the author's view, is worthy of future consideration. 

C. A FINAL NOTE 

In conclusion, this thesis has demonstrated, in effect, 
another interesting application of the method of design opti- 
mization. The author, in no way, wishes to leave the reader 
with the impression that the techniques of design optimi- 
zation are the "be all — end all" for engineering analysis. 
Design optimization techniques are useful and powerful tools 
that stand to relieve the engineer of the mundane tasks of 
numerical calculations and subsequent graphic plotting. 

But, they are just tools. In the final "analysis", good 
engineering judgment is paramount in their application and 
use . 
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APPENDIX A 



FORTRAN 


VARIABLE CROSS REFERENCE LIST 


Symbol 


Fortran Variable 


A^/Aq 


AEDVAO 




AEAOMN 


bldwt 


WEIGHT 


^.75R 


C75R 


D 


DIA 


^lim 


DIALIM 




HCL 


J 


RJ 




KQ 




KT 


noscrw 


NOSCRW 


^P 


N 




PE 




PD 


P/D 


PDIVD 


^watvap 


PWATVA 


^atm 


PATM 


promat 


PROMAT 


Qs 


QS 


^75R 


R75R 


Sc 


SC 


td 


TD 
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APPENDIX A (CONT.) 



Symbol 


Fortran Variable 


Temp 


TEMP 


(t*/c) 


TC75R 


V (ft/sec) 


V 


V (knots) 


VK 


wt 


WT 


Z 


Z 




ETAO 




ETARR 


V 


WATNU 


P 


WATRO 
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APPENDIX B 



SUBROUTINE LISTINGS 
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APPENDIX C 



ANALIZ CODES — DESIGN CASE NO. 1 
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APPENDIX E 
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APPENDIX G 



CONTROL CARD IMAGES--DESIGN CASE NO. 2 
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ANALIZ CODES — DESIGN CASE NO. 3 
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APPENDIX J 



CONTROL CARD IMAGES — DESIGN CASE NO • 3 
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